AD -A141 669 IONOSPHERIC PHYSICS(U) SACHS/FREEMAN ASSOCIATES INC 1/| 

BOWIE MD M SINGH 07 OCT 82 SFA-82-0005 
N00173-80-C-0293 


UNCLASSIFIED 


F/G 4/ 






































MICROCOPY RESOLUTION TEST CHART 

NATIONAL BUREAU OF STANDARDS-1963-A 

























PHOTOGRAPH THIS SHEET 




OD 

INVENTORY 


ert'c 'Fhyf'CS, Fm»( 

DOCUMENT IDENTIFICATION '&F4~ 22.-0005 


Noom^to-c-oiss *? 

Q I---1 


^^>TMBUT lON STATEMENT A 

A»P*OT#d fat public MboM 
Dfatribmioa Unlimited 


DISTRIBUTION STATEMENT 



DISTRIBUTION / 


AVAILABILITY CODES 


BIST AVAIL AND/OR SPECIAL 


DTIC 

k ELECTE| 

k MAY 2 111841 


DATE ACCESSIONED 


DISTRIBUTION STAMP 


84 05 17 


059 




DATE RECEIVED IN DTIC 

PHOTOGRAPH THIS SHEET AND RETURN TO DTIC-DDA-2 



DOCUMENT PROCESSING SHEET 
















AD-A141 669 


a 



FINAL REPORT 


IONOSPHERIC PHYSICS 
M. Singh, PhD 

Sachs/Freeman Associates, Inc. 
11*300 Gallant Fox Lane 
Suite 214 
Bowie, HD 20715 

07 October 1982 


SACHS/FREE MAN ASSOCIATES, INC. 

14300 GALLANT FOX LANE, SUITE 214, BOWIE, MARYLAND 20715 301-262-4400 








o 


L 



FINAL REPORT 


IONOSPHERIC PHYSICS 
M. Singh, PhD 

Sachs/Freeman Associates, Inc. 
14300 Gallant Fox Lane 
Suite 214 
Bowie, MD 20715 

07 October 1 982 


SACHS/FREEMAN ASSOCIATES, INC. 


14300 GALLANT FOX LANE, SUITE 214, BOWIE, MARYLAND 20715 301-262-4400 













L 



Security Classification: Unclassified 
SFA82-0005 


IONOSPHERIC PHYSICS 


M. Singh, PhD 

Sachs/Freeman Associates, Inc. 
14300 Gallant Fox Lane 
Suite 214 
Bowie, MD 20715 


07 October 1982 


Final Report for Period 31 July 1980 - 31 August 1982 


Distribution Statement: 
NRL Code 4187, 1230 


Prepared for: 

NAVAL RESEARCH LABORATORY 
Space Sciences Division, Code 4187 
4555 Overlook Avenue, SW 
Washington, DC 20375 


SACHS/FREEMAN ASSOCIATES, INC. 

- -: — 1 = ■ - ■ - - 1 ■ - ' — — , 


143C0 GALLANT FOX LANE, SUITE 214, 80WIE, MARYLAND 20715 301-262-4400 







REPORT DOCUMENTATION PAGE 


READ INSTRUCTIONS 
BEFORE COMPLETING FORM 


REPORT NUMBER 


SFA-RD-80-C-0293 


4 TtT\.E (end Subtitle) 


IONOSPHERIC RESEARCH 


7 AUTHORflJ 


2 GOVT ACCESSION NO. 


M. Singh, PhD 


5 Tv PE OF REPORT A PERIOD COVERED 

31 July 1980 - 31 Aug 1982 
Final Report 


6 performing org. report number 

SFA-82-0005 


a. contract or grant NUMBERS) 

N00173-80-C-0293 


12. REPORT OATE 


October 7, 1982 


I}. number of pages 


9- PERFORMING ORGANIZATION name ano address 

Sachs/Freeman Associates, Inc. 
14300 Gallant Fox Lane, Suite 214 
Bowie, Maryland 20715 


II. CONTROLLING OFFICE NAME ANO ADDRESS 

Naval Research Laboratory 

Space Sciences Div., Code 4187 

4555 Overlook Ave,, SW, Washington, DC 20375 


U. MONITORING AGENCY NAME « AOORESSCIf dlllerent from Controlling Ollice) IS. SECURITY CLASS, (ol thle report) 

DCASMA, Baltimore „ , ... , 

Room 200, 300 E. Joppa Road Unclassified 

Towson, MD 21204 


16. DISTRIBUTION STATEMENT (ot thle Report) 

Approved for public release; distribution unlimited. 


17. DISTRIBUTION STATEMENT (of the abstract entered In Block 20, If different from Report) 



19. KEY WORDS (Continue on reverse aide It necessary and Identity by block number) 

F-Region irregularities; ionospheric; S3-4 Satellite 


20 ABSTRACT (Continue on reverse side It neceseery end Identify by block number) 

The Ionospheric Diagnostic Section of the Upper Air Physics Branch, Space 
Sciences Division requires understanding similarities and differences be¬ 
tween polar and equatorial spread-F irregularities. This report discusses 
equatorial investigations with major technical, scientific and engineering 
advances contributing to the current understanding and future plans for 
adaptive communications system design and ionospheric modification and con¬ 
trol. In this report, the S3-4 satellite data analyses is summarized. 


DD i jan*73 1473 EDITION OF I NOV »S IS OBSOLETE 


Unclassifled 

SECURITY CLASSIFICATION OF THlS p'age'TwwTi 



































1.1 


INTRODUCTION 


The Upper Air Physics Branch (Code 4120) at the Naval Research Laboratory 
(NRL) conducts upper air and astrophysical research by means of rockets and 
satellites. This research includes studies of the wavelength distribution of 
solar radiation and the time variations of such emissions as well as the resonant 
scattering of solar radiation by atmospheric species. The Branch also conducts 
a rocket astronomy program to study nighttime ultraviolet and x-ray radiations 
from galactic sources and studies the effects of radiation and particles inci¬ 
dent on the Earth's atmosphere. This data is used to determine the ionic and 
molecular composition, temperature and density of the high atmosphere. The 
Ionospheric Diagnostic Section is involved with measurement of ionospheric 
plasma phenomena; specifically, F-region irregularities at equatorial and high 
latitudes. 

Accumulating information regarding the equatorial F-region ionosphere has 
proven to be important for the communications and ionospheric sciences. This 
interest has brought considerable intensity into equatorial investigations with 
major technical, scientific, and engineering advances contributing to the 
current understanding and future plans for adaptive communications system design 
and ionospheric modification and control. In this report the S3-4 satellite 
data analysis is summarized. 

2.1 EQUATORIAL IONOSPHERIC IRREGULARITIES 

The cause-effect relationship for the equatorial ionospheric irregularities 
was investigated by analyzing the S3-4 and Plumex-1 data. The results of S3-4 
satellite data analyses can be summarized as follows: 

° Near the nighttime equator there are irregularity structures with 
depletions in plasma density up to 3 orders of magnitude and range in horizontal 
extent from ten's of kms to less than 1 km. These depletions show east-west 
asymmetry with more spectral strength on the western wall of the depletions, 
suggesting more scintillations on the western wall. The mean ion mass, as infer¬ 
red from the ratio square (I /I^)^), is more inside the depletions in most of 
the depletions suggesting more molecular ions inside the depletions. 

° The power spectral analyses of S3-4 data inside the irregular structures 
show that in the intermediate wavelength domain (2km - 23m), the power spectral 
index is 2.4(±0.4). This supports the belief that these irregularities are pro¬ 
duced by collisional Rayleigh-Taylor instability. Simultaneous power spectral 
analyses on electron density and mean-ion-mass have shown a general behavior 
tending to lower power and soften spectra in ion mass fluctuations as compared 
with the fluctuations in total plasma density. 

° The ionospheric irregularities during the equatorial spread-F conditions 
are known to range from ten's (may be hundreds) of kms to 1 meter or even 11 cm. 

To investigate the plasma process(es) responsible for the generation of these 
irregularities which range six orders of magnitude in the wavelength domain, 
the histograms of spectral indices ( obtained by FFT analyses of a number of 
S3-4 orbits) in the different wavelength domains have been drawn. The long 
wavelength domain (50 km >A> 1km) lies near 1.5 while for intermediate 




(lkm >A> 20m) and transitional (100m >A> 10m) wavelengths (A> 10m) two types of 
phenomena have been observed on the power spectra: the peaks at frequencies of 
few hundreds Hz corresponding to wavelengths of few meters and a break in FFT 
near 3 meter wavelength. These are explained on the basis of excitation of 
drift wave and lower-hybrid instability respectively. 

The total spectra ranging from ten's of kms to few meters is explained by 
the concept of cascading of irregularities or heirachy of instabilities. The 
equatorial spread-F starts at wavelengths range of hundreds of kms to tens of kms 
by either TID s or gravity waves. The longer wavelengths generate the conditions 
which are ripe for the onset of Rayleigh-Taylor (R-T) instability. The gradient 
scale sizes for R-T instability is near 20 km. Once R-T instability has taken 
over, it produces the gradient scale sizes of few hundred meters which start 
gradient driven drift instability in the transitional wavelength domain. This 
phenomena of cascading of irregularities to shorter and shorter sizes produces 
conditions for excitation of drift waves, parametric instabilities and lower- 
hybrid instabilities. 

The most significant results of the investigations have been the two- 
dimensional picture of depletions and the understanding of different plasma 
processes in the different wavelength domains of ionospheric equatorial irregu¬ 
larities. 

3.1 LONGITUDINAL VARIATIONS IN EQUATORIAL IRREGULARITIES 

One of the advantages of STP/S3-4 was that it was sun-synchronous and 
polar orbit satellite which made it possible to study the longitudinal variations 
in equatorial irregularities, as well as global coverage. During the six month 
(March-September 1978) life of the satellite, about 24 percent orbits of the 
total number of orbits which covered the nighttime equator observed the ionos¬ 
pheric irregularities. Although the irregularities can be present at almost 
any longitude, 50 percent of the observed irregularities lie in the narrow 
longitude zones of 0-15° and 165-195°E. 

4.1 CALCULATION DEVELOPMENT 


To perform the above tasks, a large number of routines were developed for 
both PDP 11/10 and PDP 11/34 minicomputers. A partial list of these routines is 
given below: 

° To read the data from tapes provided by Aerojet for S3-4 data. 

° Calculation and the plot routines for the saturation currents 
(I (E), 1(1)) and the ratio square [(I3(E)/I fi (I)) 2 ] for S3-4, 
Plumex-1, and DNA/PLACES data in various formats. 

0 Plot routines for I-V characteristics for S3-4 data and the temperature 
calculation routine. 

0 Calculation of gradient scale length both for S3-4 and Plumex-1 
data and subsequent plot routines. 

0 FFT analyses routine along with the addition of window and filter. 

0 Calculation and plot routine for a chapman layer used in BIME 
simulations and SAR applications. 
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Satellite and rocket 

observations of equatorial spread-F irregularities: 
a two-dimensional model 


E. P. Szuszczewicz. M. Singh* and J. C. Holmes 
E. O Hulburt Center for Space Research. Naval Research Laboratory. Washington, DC 20375, U.S.A. 


Abstract—Recent rocket and satellite measurements of equatorial F-region irregularities have been 
able to resolve wavelengths comparable to the meter-size sensitivities of the Jicamarca and Altair 
radar backscatter techniques. In a July 1979 rocket campaign at the Kwajalein Atoll, vertical profile 
measurements by ‘in situ’ plasma probes showed the F-region marked by a number of large scale 
plasma depletions, each having its own distribution of smaller scale irregularities and a trend toward a 
co-location of the more intense irregularities with positive gradients of larger scale features. Similar 
measurements on the S3-4 Ionospheric Irregularities Satellite have shown large scale depletions (1-3 
orders of magnitude) with east-west asymmetries that point toward the western wall as the sight for 
the more intense plasma density fluctuations. The combined rocket and satellite measurements provide 
a two-dimensional model of macroscopic F-region depletions with small structures tending to develop 
more readily on the top and western boundaries. The model and associated power spectral analyses is 
in concert with a developing catalog of radar observations and the predictions of numerical simulations 
which employ the Rayleigh-Taylor instability as the primary mechanism for the generation of 
intermediate wavelength irregularities. 


1. INTRODUCTION 

Recently, there have been considerable advances in 
the understanding of equatorial spread-F as a re¬ 
sult of improved ground based radar observations 
and coordinated in situ measurements. The long 
standing Jicamarca radar results [e.g. Woodman 
and LaHoz (1976) and associated references] have 
been studied in greater detail by the expanded 
capabilities of the Altair radar facility in the Kwa¬ 
jalein Atoll (Tsunoda and Towle, 1980). These 
more recent results have investigated plume de¬ 
velopment in an east-west cross-section and pro¬ 
vided some tentative identification of plume ele¬ 
ments with local F-region plasma depletions 
(Tsunoda, 1980a, b). 

Early efforts to examine the exact relationship 
between radar plumes and ionospheric depletions 
by performing simultaneous in situ and ground- 
based radar observations (Kelley et at., 1976; 
Morse et al., 1977) were limited to conditions of 
bottomside spread-F and required extrapolations in 
space and time to establish correlations. These 
frustrations, as well as two unsuccessful attempts at 
Kwajalein in 1977 and 1978, were finally relieved 
by the successful DNA/PLUMEX campaign in 
1979 (Szuszczewicz et al., 1980a) which defini¬ 
tively established the co-location of intense radar 
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returns with the upper region of a topside F-layer 
depletion. The PLUMEX results also provided the 
first vertical profile of large-scale plasma depletions 
with superimposed distributions of smaller scale 
irregularities. 

In an attempt to establish a 2-dimensional east- 
west in situ profile that complements the develop¬ 
ing catalog of Altair observations, we present in 
this paper a composite of rocket and satellite obser¬ 
vations, from the PLUMEX campaign on the one 
hand and the S3-4 satellite on the other. These 
results are summarized in subsequent sections and 
synthesized in a way that suggests a two- 
dimensional model of macroscopic F-region deple¬ 
tions with smaller scale structures tending to de¬ 
velop more readily on the top and western bound¬ 
aries. 

2 . rocket-borne measurements 

The coordinated measurements of equatorial 
spread-F conducted during July 1979 at the Kwa¬ 
jalein Atoll involved the launch of two in¬ 
strumented rocket payloads designed to probe the 
detailed in situ structure of the turbulent ionos¬ 
pheric plasma. The first launch operation 
(PLUMEX I; 17 July 1979; 0031; 3025 LT) was 
conducted during the late phase in the development 
and decay of spread-F. The associated results and 
discussions of ground-based and rocket-borne diag¬ 
nostics have been presented by Szuszczewicz er 
al (1980a, b). 
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Fig. 1. Relative electron density profile of macroscale features as measured simultanet 
electron saturation probe currents collected on the upleg trajectory of PLUMEX I. 
intensity’ provides an approximate measure of smaller scale structure with a ±4 ir 
mately equal to ±80% fluctuations about a linear detrend (From Szuszczewicz 
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’irregularity 
•v approxi- 
■l., 1980.1 




Fig. 2. Expanded views of density fluctuations observed 


ns (\ D. H and f of Fig. 1. 
















I putposes of comparison with subst qiA m pu 
•vcl.itions of S3-4 satellite data we show here m 
I ,,. 1 the upleg measurements of relatos di istu .r 
pu-sented by correlated ion- (/„') and electron 
saturation /,‘(-V') currents collected l>\ a pair ol 
on board pulsed-plasma-probes (Hot mi s and 
Szuszczewicz. 1975: Szuszczewicz and 

HotMii-s, 19801 in the PLUMEX 1 operation. The 
ordinate has a linear scale for time-after-launch 
ssilh altitude superimposed at 50 s increments. (Be¬ 
cause ion and electron saturation currents have 
significantly different sensitivities to velocity, sheath 
and magnetic field effects, variations in l„' and 
/'(V*) not mutually corroborated were attributed 
to the various aspect sensitivities and excluded 
front Fig. 1. This approach facilitated analysis, re¬ 
duced computer time, and established credibility in 
the interpretation of the curves as relative electron 
density profiles.) 

The results in Fig. 1 show that a number of 
major depletions (ANJN, 0 ^ 0.9) were distributed 
throughout the F-region. Each of the large scale 
depletions (identified alphabetically) has its own 
distribution of irregularities, illustrated in Fig. 2 by 
the expanded view of regions C, D, H and /. It is 
clear that C is not a single narrow biteout but a 
collection of rather large irregular structures ex¬ 
tending over a total altitude domain of about 
12 km. (Vehicle velocity in region C was 2.4 km/s.) 
To develop a quantitative view of irregularity fluc¬ 
tuations observed in the F-region, contiguous 
linear detrends over 4-s intervals were executed 
throughout the entire upleg trajectory. The varia¬ 
tions about those linear detrends were then plotted 
in Fig. 1 as 'Irregularity Intensity,’ with a maximum 
relative scale of ±4. A fluctuation as great as ±4 
approximately represents a ±80% fluctuation about 
the linear detrend. (Correlation of these results 
with macroscale gradients and Altair backscatter 
contours are discussed by Szuszczewicz et al., 
1980a). The results show that the most intense 
irregularities occurred on the bottomside gradient 
(region C) with corresponding measurements at all 
other altitudes at a much lower level. We note that 
the fluctuations in the largest depletion (region 
H-l) are smaller than those at C. Furthermore, the 
fluctuations at C, D, E and I are more intense than 
at adjacent locations; and C, D, E and / are co¬ 
located with positive density gradients of large scale 
features. 

3. S3-4 SATELLITE RESULTS 
While rocket-borne instrumentation can provide 
vertical profiles of irregularities, a circular or near- 
circular orbiting satellite with high resolution in- 


xtiuiiiciil.tliiiii is required l>> assemble irregularity 
intensities ami powet law Ivlinsior in the horizontal 
plane. T he- S3-4 satellite i antes just such an experi¬ 
ment. 

The S3-4 ionospheric irregularities experiment 
employs a pair of pulsed-plasma-probes (Holmes 
and Szuszczewicz, 1975; Szuszczewicz and 
Holmes, 1980) on a polar (96.4° inclination), sun- 
synchronous (2230 h LT. equatorial crossing), F- 
region orbiting satellite. The experiment provides 
direct measurements of the ionospheric state 
{N„ T,), its condition of irregularity (6N,), and 
associated electron density fluctuation power 
spectra [P(fc)] with 5-20 m resolution. 

Figure 3 presents a sample of relative electron 
density data collected during a nighttime equatorial 
crossing with the satellite orbiting at an altitude of 
240 km with a velocity of 7.53 km/s. The velocity 
component pcrpendict 'ar to the geomagnetic field 
is 2.0 km/s. 50 s of data are displayed covering a 
total horizontal extent equal to 366 km, moving 
from east-to-west with increasing time. 

While the data represents a transit from the 
southeast to the northwest, our discussions and 
associated interpretations will assume that all varia¬ 
tions along the magnetic field are negligible. We 
therefore interpret our results in the framework of 
a two-dimensional system, perpendicular to the 
magnetic field. This assumption is similarly em¬ 
ployed in theoretical analyses of equatorial spread- 
F (e.g. Scannapieco and Ossakow 1976) and 
supported by 6300 A airglow intensity measure¬ 
ments of F-region depletions which show signific¬ 
antly greater variation in east-west extent than 
north-south (Weber et al., 1978). Our discussions 
will also assume temporal variations to be nonexis¬ 
tent during the satellite traversal across any given 
depletion (t'-pically 10 s or less in Fig. 3). We 
consider this a reasonable position since saturation 
times for unstable Rayleigh-Taylor modes are of 
the order of 4000 s (Keskinen et al, 1980a) and 
F-layer vertical drift velocities are generally in the 
10-20 m/s range (e.g. Szuszczewicz et al., 1980). 
In neither case should our interpretations be in¬ 
fluenced by a 10 s time lapse as the satellite passes 
from the eastern to the western boundary of a 
depletion. 

The data reveal four large scale depletions rang¬ 
ing up to 98 km wide and a factor of 300 in depth, 
with smaller scale structure visible down to a frac¬ 
tion of a kilometer. The upper portion of Fig. 3 
presents 'irregularity intensity’ as measured by rela¬ 
tive fluctuations about continuous 0.27 s (2.1km) 
linear detrends. 

Further discussion is facilitated by identifying 
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Fig 3. The electron density variations as measured by the electron probe aboard S3-4 on rev No 
2123 On the top is the irregularity intensity (81/1) determined over contiguous 2.1 km intervals 
throughout the depletions 


certain features in Fig. 3. First, there are clearly 
'defined regions of undisturbed background ionos¬ 
phere, marked alphabetically A through D; the 
smoothness of the relative density and the corres¬ 
ponding 0% fluctuations attest to their undisturbed 
nature. The eastern boundary is defined as the 
region of density gradient moving westward from 
an undisturbed domain to the least lower bound of 
relative plasma density in the depletion. Everything 
to the west of that minimum is defined as the 
western boundary. Admittedly these definitions 
leave no room for a region that might be called the 
‘depletion center,' but depletions 3 and 4 suggest 
that there is no easily defined ‘depletion center.’ 
This is also borne out by other sets of P 3 data. 

Focussing on depletions 3 and 4, we see that the 
irregularity intensities are 2 to 3 times larger on the 
western boundary than on its eastern counterpart. 
This same relationship is true in depletions 1 and 2, 
but only after a qualification that suggests that 1 
and 2 are halves of a larger depletion bounded by 
A and B. This is supported in part by the nonexis¬ 
tence of a quiescient ionosphere between the two. 
When viewed from this perspective the western 
boundary is approximately twice as intense in ir¬ 
regularity intensity as the eastern boundary. 

Signatures in irregularity strengths and relation¬ 
ships to plasma instability mechanisms c ut he 


further explored through power spectral density 
analyses. We present in Fig. 4 just such results for 
each of the boundaries in Fig. 3 (IE and 1W refer 
respectively to the eastern and western boundary of 
depletion number 1). Power spectral analyses are 
presented across the boundaries of each of the four 
depletions with spectral indices (q, in a spectral fit 
to P = P 0 f"') ranging from 1.9 to 2.3. (These indi 
ces are in agreement with the work of Kfskinkn et 
al. (1980) where numerical simulations predict n 
2.0 to 2.5 for horizontal irregularity structures per¬ 
pendicular to B. This result runs parallel to an 
earlier comparison (Sec s/t zt witz and Hoi.su s. 
1980; Kesmnen el al. 1980b) which showed that 
the power spectral index across region C in 
PLUMEX 1 also agreed with the predictions of 
Keskinen el al. (1980a) in the vertical plane 1 
More important to the discussion of east-west 
asymmetries is the ratio of the spectral strength' 
(P„) across the two boundaries. By defining Pi as 
the P 0 (West)/P 0 (East), we find the ratio P, to ex¬ 
tend from 1.4 to 11.3, i.e. the irregularity spectral 
strength on the western wall is 1.4 to 11.3 times 
more intense than on the east.ern counterpart. 

These results are in concert with the radar obser¬ 
vations of Tsunooa (1979) which show that the 
bottomside backscatter strength is often asymmet* 
ric in the east-west plane, with the western side of 
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Fig. 4. Power spectral analyses across the boundaries of the four depictions (1-4) in Fig. 3. The 
number in the upper portion of each | <el is the spectral index N. in a P = P t J " fit to the data. The 
value P, shown in each right hand pan orresponds to the ratio of spectral strengths P, ■* P 0 (western 
boundary) - P„ (eastern boundary). 


the plumes being the more intense. The combined 
observations support a model of E-W asymmetry 
which allows for a neutral-wind-driven instability- 
growth-rate enhancement on the western side of a 
rising bottomside F-region depletion (Tsunoda, 
1979). 

4. CONCLUSIONS 

The combination of recent rocket and satellite 
data allows for the development of a two- 
dimensional empirical model of equatorial deple¬ 
tions and associated irregularities. Elements in the 
data have led to the following conclusions: 


(1) Equatorial depletions are macrostructures 
within which much smaller scale irregularities are 
imbedded; 

(2) These smaller scale irregularities (in the ver¬ 
tical plane) tend to derive their energies from posi¬ 
tive density gradients on the topside of local deple¬ 
tions; 

(3) Vertical and horizontal power spectral anal¬ 
yses display indices from 1.8 to 2.5, a result which 
is consistent with Rayleigh-Taylor instability gener¬ 
ation of intermediate wavelength irregularities dur¬ 
ing the occurrence of equatorial spread-F (Keskj- 
nen er al. 1980a); 
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(4) When viewed horizontally, bottomside deple¬ 
tions have an east-west asymmetry with the more 
intense fluctuations and spectral strengths being 
observed on the western boundary. This result is 
consistent with the radar measurements of 
Tsunoda (1979) and scintillation observations of 
Livingston et al. (1980): 


(5) I lie combination of all results leads to a 
two-di -nsional model of macroscopic T-reeion 
depletions with smaller scale structures developing 
on the top and western boundaries. This model and 
associated spectral indices fully supports the com¬ 
putational work of Keskinen et al. (1980a). 
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associated with a density gradient and/or current-driven plasma instability. We have 
made preliminary analyses of the density fluctuations and FFT power spectra for 
evidence of characteristic scale sizes and power law dependence. Scale sizes from 
10-300 km are clearly evident in the irregularities; power law fits to the spectra have 
spectral indices in the range -1.5 to -1.7. For the diffuse aurora our results suggest 
support for a recent numerical study of the nonlinear evolution of the current 
convective instability. 
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THE STP/S3-4 SATELLITE EXPERIMENT: 

HIGH LATITUDE LARGE SCALE DENSITY IRREGULARITIES 


INTRODUCTION 

The polar regions of the earth’s ionosphere are known 
to exhibit a wide range of plasma kinetic interactions that 
are very dynamic and complex. The plasma kinetic processes 
have interaction scale sizes that range from the order of 
auroral longitudinal dimensions down to particle gyroradii. 
Although the phenomenology and morphology of auroral processes 
are well known, the physical mechanisms are not well under¬ 
stood. Plasma density irregularities in the auroral F- 
region have been detected through ground-based scintillation 
observations (Aarons, 1973; Fremouw et al., 1977; Buchau et 
al., 1978; Rino et al., 1978), satellite measurements (McClure 
and Hanson, 1973; Dyson et al., 1974; Sagalyn et al., 1974; 
Phelps and Sagalyn, 1976; Clark and Raitt, 1976; Weber and 
Buchau, 1981), and rocket observations (Kelley et al., 

1980). Several physical mechanisms have been suggested to 
explain the high latitude irregularities. A recent review 
(Fejer and Kelley, 1980) lists three general sources of the 
irregularities: 1) particle production (especially by low 
energy precipitating electrons), 2) electrostatic turbulence, 
and 3) plasma instabilities. Except in extremely simple 
situations, it is likely that all three mechanisms will 
contribute to the spectrum of irregularities. 

In this report we will examine two cases of high latitude 
measurements of electron density irregularities obtained in 
the nightside auroral F-region with the pulsed plasma probe 
experiment on the STP/S3-4 satellite. These two cases 
suggest the importance of plasma instabilities associated 
with particle precipitation, i.e, field-aligned currents, 
and illustrate the variety of density irregularities en¬ 
countered at high latitudes. The first observation we 
discuss is that associated with a quiescent aurora, in which 
the structure of density irregularities is relatively well 
defined on several scale lengths, from 10 km to about 300 
km. These scale lengths are significant for some recent 
theoretical studies on the applicability of the current 
convective instability in diffuse auroras. The second 
observation we discuss was obtained during a more dynamic 
state of the aurora and serves to illustrate some of the 
complexity associated with possible gradient- and/or current- 
driven plasma instabilities in discrete auroras. 

MEASUREMENTS AND ANALYSIS 

Plasma Density Measurements 

The polar-orbiting STP/S3-4 satellite is instrumented 
with a pair of pulsed-plasma-probes that allow high resolution 
measurements of electron density N , temperature T g , density 
fluctuation power spectra P (k), plasma potential 6 V and 
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variations in the mean ion mass 6M./M. (Szuszczewicz et al . , 
1981]. We briefly review the experimental operation only 
for the particular data mode of interest here. In the low- 
data-rate mode, the double probe system sampled electron- 
and ion- saturation currents simultaneously at 100 samples 
s 1 , Under conditions of fixed-payload potential and stable 
velocity-field and magnetic-field aspect (see e.g., Szuszczewicz 
and Takacs, 1979) , the density N and density fluctuation 
AN are directly proportional to e the probe saturation current 
I e and its fluctuations AI . The STP/S3-4 satellite was in 
a e sun synchronous orbit ^ith an inclination of 96.5° 
and crossed the equatorial plane at local times of 2230 and 
1030. The altitude ranged from 160 to 260 km. In the 
discussion below of irregularity scale lengths, we will use 
an average spacecraft velocity of 8 km s 1 to convert time 
variations to spatial variations, using AL=vAt, where AL is 
the spatial length, v is the spacecraft velocity, and At is 
the sampled time interval. 

Quiescent Aurora 

About five minutes of saturation currents, I and I., 
obtained on orbit 390 are shown in Figure 1. These measurements 
were obtained in the night side auroral F-region over the 
south pole at an altitude of 260 km with the spacecraft 
moving from the south polar cap region toward the nightside 
equator (left to right in Figure 1). The probe currents 
clearly show that irregularities were detected on both 
electron and ion probes. The detailed tracking of one probe 
by the other indicates that the probe system was working 
properly and that the irregularities were due primarily to 
variations in the ambient electron and ion densities and not 
the result of secondary effects such as changing spacecraft 
potential, aspect sensitivity, etc. In this preliminary 
report, only the electron density irregularities will be 
discussed, and we simply note the correlative behavior of 
the ion density irregularities. 

The large depletion in electron density (almost an 
order of magnitude) centered at 1840 UT corresponds to a 
diffuse aurora. The geomagnetic Kp index for the period of 
this measurement was 1+, indicating relatively quiescent 
auroral conditions. The ion density shows a similar depletion. 
Both the electron and ion depletions extend over a time 
interval of about 25 s, corresponding to a distance of about 
200 km along the approximately south-north direction of 
spacecraft motion. This aurora was also detected by Huffman 
et al. (1980) with the vacuum ultraviolet (VUV) experiment 
on board the same spacecraft. A detailed comparison of the 
density measurements in Figure 1 with the VUV measurements 
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Fig. 1 - The electron and ion currents sampled in the nightside 
auroral F-region during a pass through a quiescent diffuse aurora 









(Fig. 15 in Huffman et al.) indicates that a good correlation 
exists. A few discrete arcs detected by the VUV measurement 
in the time interval 1839-1840 UT appear to be correlated « 

with some of the density structure. Equatorward (to the 
right) of the large depletion, the electron density shows a 
relative enhancement with a scale length of about 300 km. 

A more detailed plot of the electron density depletion 
of Figure 1 and its contiguous domains is shown in Figure 2. 

(The expanded interval is indicated by the arrows above the 
electron current trace in Figure 1.) In this plot every 
data point sampled at 100 s 1 is shown and the density scale 
on the left of the figure is approximate. Within the depletion 
at least three well-defined density variations with nearly 
periodic behavior are indicated. Assuming the auroral structure 
to be relatively stable, the width of the three density 
structures along the spacecraft trajectory is about 20 km 
and the large scale depletion has a width of about 200 km. 

On either side of the large scale depletion the density 
irregularities generally appear to have high and low frequency 
components with the low frequency components having scale 
widths comparable to the widths of the density structures 
within the depletion, i.e., about 20 km. 

The Fast Fourier Transform (FFT) power spectrum of the 
density irregularities in each of three intervals are shown 
in the insets to Figure 2. The data time interval over 
which the FFT was calculated is indicated below each spectrum. 

The Nyquist frequency for each spectrum is 6.25 Hz, corresponding 
to a minimum scale size of about 0.6 km. Each spectrum is 
normalized to the maximum power calculated from the FFT. 

The straight line drawn through each spectrum is the least- 
squares power law P=P f , where P/P is the normalized 
power, f is the frequency, and n is°the spectral index. 

The spectra are calculated without smoothing, except for 
the removal of the lowest frequency (DC) component by applying 
a linear detrend. 

The first power spectrum (on the left) suggests a k- 
dependence of k 1 ' “ 8 , where k = 2 tt/X is the wave number and 
X is the wavelength of the spatial irregularity. The value 
of the index is consistent with previous reports (Dyson et 
al., 1974; Sagalyn et al., 1974; Phelps and Sagalyn, 1976; 

Kelley et al., 1980) on high latitude irregularities where 
scale sizes from less than 1 km up to 200 km have been 
considered. In the second power spectrum, taken within the 
large scale depletion, the spectral index implies a k 1,58 
Variation in spatial irregularity. The low frequency 
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S3-4 ORBIT 390 NIGHTSIDE AURORAL REGION 9 APRIL 1978 



UT IHRMN : SC) 1839 : 30 1840 : 00 1840 : 30 

MLAT (D*g) -64.7° -62.8° -80.7° 

Fig. 2 - An expanded time-scale plot of electron density irregularities in 
the diffuse aurora of Figure 1. The power spectra indicate the frequency 
dependence of the least-squares power law. Frequency f is directly propor¬ 
tional to wave-number k, assuming the irregularities are spatial structures. 
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spectral range 0.1-0.3 Hz corresponds to the well-defined 
structures of about 20 km half-width in the density-vs- 
time plot. The spectral power in the 20-km structures 
inside the large density depletion is about three times 
greater than the power in similar scale size structures 
outside the depletion, a fact not obvious in the power 
spectra because of the normalization. The third power 
spectrum shows a power law dependence with a spectral index 
of -1.69. 

An examination of the density-vs-time plot of Figure 2 
shows that approximately 20-km scale size structures are 
clearly evident both inside and outside the large depletion. 

A possible source of these 20-km density irregularities is 
suggested by recent theoretical studies on the current 
convective instability in diffuse auroras (Ossakow and 
Chaturvedi, 1979; Keskinen et al., 1980; Huba and Ossakow, 
1980). The current convective instability is proposed as a 
mechanism for generating large scale ('v 10-50 km) density 
irregularities in the diffuse auroral F-region. The current 
convective instability may thus provide an explanation for 
the sheet-like density irregularities that have been detected 
by Fremouw et al. (1977) and Rino et al. (1978) in their 
analysis of enhanced scintillations at auroral latitudes. In 
a numerical analysis, Keskinen et al., (1980) considered the 
nonlinear evolution of the current convective instability 
with an initial density gradient scale size of 50 km. They 
obtained a spectrum of irregularities with a spectral index 
of -1.7 in the north-south direction. This scale size and 
spectral index are similar to the scale lengths and to the 
spectral indices shown in Figure 2. Thus we suggest that 
the observations provide experimental support for the current 
convective instability as a possible mechanism for the 
density irregularities associated with the diffuse aurora. 

Dynamic Aurora 

The plasma probe measurements of Figure 3 illustrate 
the case of a dynamic aurora observed on orbit 244 in the 
southern night-side auroral F-region. The K index for this 
period was 3+, indicating that moderately active geomagnetic 
perturbations were occurring. As in Figure 1, the electron 
and ion saturation currents were sampled in the low- 
data-rate mode (100 s *)• The two currents track each other 
very well so that we are confident of proper measurements 
directly proportional to the ambient plasma density. The 
outstanding features of the density measurements in Figure 3 
are the sharp transition in density occurring at about 
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Fig. 3 - Hie electron and ion currents sampled in the nightside 
auroral F-region for a dynamic state of the aurora. 


7 







1915:15 UT and the large peak-to-peak variations between 
1916 and 1917 UT. The sharp transition in density is more 
apparent in the electron current than in the ion current 
because of noise interference occurring in the ion probe 
before about 1915:10 UT. A detailed comparison of the two 
currents shows that the transition itself is unaffected by 
the noise and the two currents track each other almost 
exactly within the transition. 

A detailed plot of the electron current in the time 
interval indicated by the arrows in Figure 3 is shown in 
Figure 4 with an approximate density scale to the left. The 
sharp transition begins at about 1915:17 UT and consists of 
large irregular oscillations over the transition width of 
about 1.5 s. The density jump associated with the transition 
is about a factor of 4. Immediately behind the transition 
there are at least four distinct oscillations with an average 
period of about 2 s. These oscillations seem to decay 
exponentially, as indicated by the dashed line (on a log 
scale, an exponential curve will appear as a straight line). 

The VUV observations on the same spacecraft indicate that 
auroral emissions were detected at the same time the density 
transition occurred (R.E. Huffman, private communication) 
suggesting that the density jump is associated with a discrete 
auroral arc. The observation of a sharp transition suggests 
that a moving boundary has been detected and we conclude 
that the density jump is associated with the surge of an 
auroral arc. The surge direction cannot actually be determined 
from our observations; generally the movement is westward 
and poleward (Akasofu, 1968). The sharp transition in 
density implies that gradient-driven plasma instabilities 
are likely to occur; with our interpretation of a surging 
auroral arc, the decaying oscillations behind the density 
jump would represent spatial oscillations left behind by the 
movement of the source of instability, the surge front. The 
electron density increase by a factor of 4 across the surge 
front probably results from the enhanced ionization due to 
precipitating particles. Similar ionization enhancements of 
the F-layer associated with discrete auroral arcs have been 
observed with ionospheric sounders (Pike et al., 1977). The 
surge front very likely has field-aligned currents associated 
with it (Iijima and Potemra, 1978), so that current-driven 
plasma instabilities are also likely to be present and may 
contribute to the density oscillations. It is possible that 
the large and irregular variations seen at later times 
(after 1915:55 UT) also may have had their origin in the 
surge of the auroral arc. If the irregular density structures 
behind the surge front are quasi-stationary, then a wide 
range of scale lengths is evident, from about 4 km to 70 km 
and longer. In contrast, the region ahead of the surge 








ELECTRON DENSITY (cm- ; 


S3-4 ORBIT 244 NIGHTSIDE AURORAL REGION 31 MARCH 1978 



UT (HRMN : SC) 1915 : 00 1915 : 30 1916 : 00 

MLAT(Dbq) -75.1° -73.2° -71.3° 

Fig. 4 - An expanded time-scale plot of the electron density for the 
aurora of Figure 3. The sharp density jump is identified as the surge 
of a discrete auroral arc. An exponential decay for the well-defined 
oscillations behind the surge front is suggested by the dashed line. 
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front has a distinctly different appearance, with only low 
amplitude, high frequency variations. This region ahead of 
the surge front may be a relatively quiescent ionosphere 
about to be perturbed by the surging auroral arc. The 
different character of the irregularities on either side of 
the density jump that we have identified as a surging 
auroral arc is additional evidence that the density jump is 
a dynamic transition between two states of the F-region 
ionosphere. 

The power spectra for three time intervals are also 
shown in Figure 4. Just ahead of the density jump the 
spectrum of the irregularities shows that the slope of the 
least-squares power law is relatively flat, suggesting that 
thermal fluctuations dominate the spectrum and the ionosphere 
is in a quiescent state. The spectrum immediately behind 
the density jump shows a power law with steeper slope, 
indicating the presence of enhanced power at low frequencies. 
The well-defined decaying oscillations behind the density 
jump appear in the spectrum as the component centered on a 
frequency of about 0.3 Hz. Under the assumption that they 
are spatial variations, these oscillations have a scale 
length of about 10 km. The spectral index n = -1.66 for the 
second spectrum in Figure 4 is nearly identical to one of 
the spectra in Figure 2. This near-equality suggests that 
similar transfer processes in spectral energy may be involved. 
(We note, however, that a specific spectral index does not 
uniquely define an instability process.) The third spectrum 
in Figure 4 corresponds to the. region of large peak-to-peak 
variations. The spectrum clearly shows the predominance of 
structure in the range 0.1-1.0 Hz (40-4 km). The spectral 
index is near -2, possibly due to the presence of steep 
edges in the irregularities. 

CONCLUDING COMMENTS 


The aurora-related F-region density irregularities 
discussed in this report illustrate the variety of plasma 
kinetic processes that are involved in high latitude ionospheric 
effects. Our study of the density irregularities is in a 
preliminary stage, and we expect that new and significant 
results will be forthcoming as the study advances. It is 
already obvious from an initial survey of available satellite 
data that aurora-related effects are detected outside the 
average auroral boundaries, i.e., both at higher and lower 
latitudes. Since the high latitude region is directly 
associated with magnetosphere-ionosphere coupling processes, 
we anticipate that the analysis of the ionospheric density 
irregularities in regions such as the polar cusp will address 
problems of global importance. 
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FT. HUACHUCA, AZ 85613 
01CY ATTN CCC-EMEO GEORGE LANE 

WWMCCS SYSTEM ENGINEERING ORG. 
WASHINGTON, D.C. 20305 
01CY ATTN R. CRAWFORD 

COMMANDER/DIRECTOR 
ATMOSPHERIC SCIENCES LABORATORY 
U.S. ARMY ELECTRONICS COMMAND 
WHITE SANOS MISSILE RANGE, NM 88002 
01CY ATTN OELAS-EO F. NILES 

OIRECTOR 

BMD ADVANCED TECH CTR 
HUNTSVILLE OFFICE 
P.O. BOX 1500 
HUNTSVILLE, AL 35807 
01CY ATTN ATC-T MELVIN T. CAPPS 

01CY ATTN ATC-0 W. DAVIES 

01CY ATTN ATC-R DON RUSS 

PROGRAM MANAGER 
BMD PROGRAM OFFICE 
5001 EISENHOWER AVENUE 
ALEXANDRIA, VA 22333 
01CY ATTN DACS-BMT J. SHEA 

CHIEF C-E SERVICE DIVISION 
U.S. ARMY COMMUNICATIONS CMD 
PENTAGON RM. 1B269 
WASHINGTON, D.C. 20310 
01CY ATTN C-E-SERVICES DIVISION 

COMMANDER 

FRADCOM TECHNICAL SUPPORT ACTIVITY 
DEPARTMENT OF THE ARMY 
FORT MONMOUTH, N.J. 07703 
01CY ATTN ORSEL-NL-RO H. BENNET 
01CY ATTN ORSEL-PL-ENV H. BOMKE 
01CY ATTN J.E. QUIGELY 


COMMANDER 

HARRY DIAMOND LABORATORIES 
DEPARTMENT OF THE ARMY 
2800 POWDER MILL ROAD 
ADELPHI, MD 20783 

(CNWOI-INNER ENVELOPE: ATTN: DELHD-RBH) 
01CY ATTN DELHD-TI M. WEINER 

01CY ATTN DELHO-RB R. WILLIAMS 

01CY ATTN DELHD-NP F. WIMENITZ 

01CY ATTN DELHD-NP C. MOAZED 

COMMANDER 

U.S. ARMY COMM-ELEC ENGRG INSTAL AGY 
FT. HUACHUCA, AZ 85613 
01CY ATTN CCC-EMEO GEORGE LANE 

COMMANDER 

U.S. ARMY FOREIGN SCIENCE & TECH CTR 
220 7TH STREET, NE 
CHARLOTTESVILLE, VA 22901 
01CY ATTN DRXST-SD 

01CY ATTN R. JONES 

COMMANDER 

U.S. ARMY MATERIAL DEV & READINESS CMD 
5001 EISENHOWER AVENUE 
ALEXANDRIA, VA 22333 
01CY ATTN DRCLDC J.A. BENDER 

COMMANDER 

U.S. ARMY NUCLEAR AND CHEMICAL AGENCY 
7500 BACKLICK ROAD 
BLDG. 2073 

SPRINGFIELD, VA 22150 
01CY ATTN LIBRARY 

DIRECTOR 

U.S. ARMY BALLISTIC RESEARCH LABS 
ABERDEEN PROVING GROUND, MD 21005 
01CY ATTN TECH LIB EDWARD BAICY 

COMMANDER 

U.S. ARMY SATCOM AGENCY 
FT. MONMOUTH, NJ 07703 
01CY ATTN DOCUMENT CONTROL 

COMMANDER 

U.S. ARMY MISSILE INTELLIGENCE AGENCY 
REDSTONE ARSENAL, AL 35809 
01CY ATTN JIM GAMBLE 




DIRECTOR 

U.S. ARMY TRAOOC SYSTEMS ANALYSIS ACTIVITY 
WHITE SANDS MISSILE RANGE, NM 88002 
01CY ATTN ATAA-SA 
01CY ATTN TCC/F. PAYAN JR. 

01CY ATTN ATAA-TAC LTC. J. HESSE 

COMMANDER 

NAVAL ELECTRONIC SYSTEMS COMMAND 
WASHINGTON. D.C. 20360 
01CY ATTN NAVLEX 034 T. HUGHES 
01CY ATTN PME 117 
01CY ATTN PME 117-T 
01CY ATTN CODE .5011 
01CY ATTN PME-106-T 

COMMANDING OFFICER 
NAVAL INTELLIGENCE SUPPORT CTR 
4301 SUITLAND ROAD, BLDG. 5 
WASHINGTON, O.C. 20390 
OlCY, ATTN MR. DUBBIN STIC 12 
OlCY ATTN NISC-50 
OlCY ATTN CODE 5404 J. GALET 

COMMANDER ! 

NAVAL SURFACE WEAPONS CENTER 
DAHLGREN LABORATORY 
DAHLGREN, VA 22448 
OlCY ATTN CODE DF-14 R. BUTLER 

OFFICE OF NAVAL RESEARCH 
ARLINGTON, VA 22217 
OlCY ATTN CODE 465 
OlCY ATTN CODE 461 
OlCY ATTN COOE 420 
OlCY ATTN CODE 421 

COMMANDER 

AEROSPACE DEFENSE COMMAND/DC 
DEPARTMENT OF THE AIR FORCE 
ENT AFB, CO 80912 
OlCY ATTN DC MR. LONG 

COMMANDER 

AEROSPACE DEFENSE COMMAND/XPD 
DEPARTMENT OF THE AIR FORCE 
ENT AFB, CO 80912 
OlCY ATTN XPDQQ 
OlCY ATTN XP 


AIR FORCE Cl Of’HYSICS LABORATORY 
HANSCOM APB, MA 01731 
OlCY ATTN OPR HAROLD GARDNER 
OlCY ATTN OPR-1 JAMES C. ULWICK 
OlCY ATTN LKB KENNETH S. W. CAMPION 
OlCY ATTN OPR ALVA T. STAIR 
OlCY ATTN PHD JURGEN BUCHAU 
OlCY ATTN PHD JOHN P. MULLEN 

AF WEAPONS LABORATORY 
KIRTLAND AFB, NM 87117 
OlCY ATTN SUL 

OlCY ATTN CA AUTHER H. GUENTHER 
OlCY ATTN NYTC CAPT. J. BARRY 
OlCY ATTN NYTC JOHN M. KAMM 
OlCY ATTN NYTC CAPT MARK A. FRY 
OlCY ATTN NTES MAJ GARY GANONG 
OlCY ATTN NYTC J. JANNI 

AFTAX 

PATRICK AFB, FL 32925 
OlCY ATTN TF/MAJ WILEY 
OlCY ATTN TN 

AIR FORCE WRIGHT AERONAUTICAL LABS 
WRIGHT-PATTERSON AFB, OH 45433 
OlCY ATTN AMD, WADE HUNT 
OlCY ATTN AMD, ALLEN JOHNSON 

DEPUTY CHIEF OF STAFF 
RESEARCH, DEVELOPMENT, S ACQ 
DEPARTMENT OF THE AIR FORCE 
WASHINGTON, D.C. 20330 
OlCY ATTN AFRDQ 

HEADQUARTERS 

ELECTRONIC SYSTEMS DIVISION/XR 
DEPARTMENT OF THE AIR FORCE 
HANSCOM AFG, MA 01731 
OlCY ATTN XR J. DEAS 

HEADQUARTERS 

ELECTRONIC‘SYSTEMS D1VISION/YSEA 
DEPARTMENT OF THE AIR FORCE 
HANSCOM AFB, MA 01731 
OlCY ATTN YSEA 

COMMANDER 

NAVAL OCEAN SYSTEMS CENTER 
SAN DIEGO, CA 92152 
03CY ATTN CODE 532 W. MOLER 
OlCY ATTN CODE 0230 C. BAGGETT 
OlCY ATTN CODE 81 R. EASTMAN 
OlCY ATTN CODE 2200 H. RICHTER 
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OIRECTOR 

NAVAL RESEARCH LABORATORY 
WASHINGTON. D.C. 20375 
01CY CODE 4100 
01CY CODE 4101 
01CY CODE 4120 
01CY CODE 4701 JACK D. BROWN 
20CY CODE 2628 
01CY CODE 4732 E. MCLEAN 
01CY CODE 6000 
01CY CODE 7000 
01CY CODE 7500 
01CY COOE 7580 
01CY CODE 7551 
01CY COOE 7555 
01CY CODE 7900 

COMMANDER 

NAVAL SEA SYSTEMS COMMAND 
WASHINGTON, D.C. 20362 
01CY ATTN CAPT. R. PITKIN 

COMMANDER 

NAVAL SPACE SURVEILLANCE SYSTEM 
DAHLGREN, VA 22448 
01CY ATTN CAPT. J.H. BURTON 

OFFICER-IN-CHARGE 
NAL SURFACE WEAPONS CENTER 
WHITE OAK, SILVER SPRING, MD 20910 
01CY ATTN CODE F31 

DIRECTOR STRATEGIC SYSTEMS PROJECT OFFICE 
DEPARTMENT OF THE NAVY 
WASHINGTON, D.C. 20376 
01CY ATTN NSP-2141 
01CY ATTN NSSP-2722 FRED WIMBERLY 

NAVAL SPACE SYSTEM ACTIVITY 
P.O. BOX 92960 
WORLDWAY POSTAL CENTER 
LOS ANGELES, CA 90009 
01CY ATTN LCOR DONALD SNOOOY 
01CY ATTN COMMANDING OFFICER 

HEADQUARTERS 

ELECTRONIC SYSTEMS DIVISION/DC 
DEPARTMENT OF THE AIR FORCE 
HANSCOM AFB, MA 01731 
01CY ATTN OCKC MAJ J.C. CLARK 


COMMANDER 

FOREIGN TECHNOLOGY DIVISION, AFSC 
WRIGHT-PATTERSON AFB, OH 45433 
01CY ATTN NICD LIBRARY 
01CY ATTN ETOP B. BALLARD 

COMMANDER 

ROME AIR DEVELOPMENT CENTER, AFSC 
GRIFFISS AFB, NY 13441 
01CY ATTN DOC LIBRARY/TSLD 
01CY ATTN OCSE V. COYNE 

SAMSO/SZ 

POST OFFICE BOX 92960 
WORLDWAY POSTAL CENTER 
LOS ANGELES, CA 90009 
(SPACE DEFENSE SYSTEMS) 

01CY ATTN SZJ 

STRATEGIC AIR COWAND/XPFS 
OFFUTT AFB, NB 68113 
01CY ATTN XPFS MAJ. B. STEPHAN 
01CY ATTN ADWATE MAJ. BRUCE BAUER 
01CY ATTN NRT 

01CY ATTN DOK CHIEF SCIENTIST 

SAMSO/YA. 

P.O. BOX 92960 
WORLDWAY POSTAL CENTER 
LOS ANGELES, CA 90009 
01CY ATTN YAT CAPT L. BLACKWELDER 

SAMSO/SK 
P.O. BOX 92960 
WORLDWAY POSTAL CENTER 
LOS ANGELES, CA 90009 

01CY ATTN SKA (SPACE COMO SYSTEMS) 

M. CLAVIN 

SAMSO/MN 

NORTON AFB, CA 92409 
(MINUTEMAN) 

01CY ATJN MNNL LTC KENNGDY 
COMMANDER 

ROME AIR DEVELOPMENT CENTER, AFSC 
HANSCOM AFB, MA 01731 

01CY ATTN EEP A. LORENTZEN 
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DEPARTMENT OF ENERGY 

DEPARTMENT OF ENERGY 
ALBUQUERQUE OPERATIONS OFFICE 
P.O. BOX 5400 
ALBUQUERQUE, NM 87115 
01CY ATTN OOC CON FOR 0. SHERWOOD 

DEPARTMENT OF ENERGY 
LIBRARY ROOM G-042 
WASHINGTON, O.C. 20545 
01CY ATTN DOC CON FOR A. LABOWITZ 

EG&G, INC. 

LOS ALAMOS DIVISION 
P.O. BOX 809 
LOS ALAMOS, NM 85544 
01CY ATTN DOC CON FOR J. BREEDLOVE 

UNIVERSITY OF CALIFORNIA 
LAWRENCE LIVERMORE LABORATORY 
P.O. BOX 808 
LIVERMORE, CA 94550 
01CY ATTN DOC CON FOR TECH INFO DEPT 

01CY ATTN DOC CON FOR 1-389 R. OTT 

01CY ATTN DOC CON FOR L-31 R. HAGER 

01CY ATTN DOC CON FOR L-46 F. SEWARD 

LOS ALAMOS SCIENTIFIC LABORATORY 
P.O. BOX 1663 
LOS ALAMOS, NM 87545 
01CY ATTN DOC CON FOR J. WOLCOTT 

01CY ATTN DOC CON FOR R.F. TASCHEK 

01CY ATTN DOC CON FOR E. JONES 

01CY ATTN DOC CON FOR J. MALIK 

01CY ATTN DOC CON FOR R. JEFFRIES 

01CY ATTN DOC CON FOR J. ZINN 

01CY ATTN DOC CON FOR P. KEATON 

01CY ATTN DOC CON FOR D. WESTERVELT 

01CY' ATTN DOC CON FOR M. PONGRATZ 

SANDIA LABORATORIES 
P.O. BOX 5800 
ALBUQUERQUE, NM 87115 
01CY ATTN DOC CON FOR J. MARTIN 

01CY ATTN DOC CON FOR W. BROWN 

01CY ATTN DOC CON FOR A. THORNBROUGH 

01CY ATTN OOC CON FOR T. WRIGHT 

01CY ATTN OOC CON FOR D. DAHLGREN 

01CY ATTN DOC CON FOR 3141 

01CY ATTN DOC CON FOR SPACE PROJECT OIV 


SANDIA LABORATORIES 
LIVERMORE LABORATORY 
P.O. BOX 969 
LIVERMORE, CA 94550 
01CY ATTN DOC CON FOR B. MURPHY 
01CY ATTN DOC CON FOR T. COOK 

OFFICE OF MILITARY APPLICATION 
DEPARTMENT OF ENERGY 
WASHINGTON, D.C. 20545 
01CY ATTN DOC CON FOR D. GALE 

OTHER GOVERNMENT 

CENTRAL INTELLIGENCE AGENCY 
ATTN RD/SI, RM 5G43, HQ BLDG 
WASHINGTON, D.C. 20505 
01CY ATTN OSI/PSID RM 5F 19 

DEPARTMENT OF COMMERCE 
NATIONAL BUREAU OF STANDARDS 
WASHINGTON, O.C. 20234 

(ALL CORRES: ATTN SEC OFFICER FOR) 

01CY ATTN R. MOORE 

DEPARTMENT OF TRANSPORTATION 
OFFICE OF THE SECRETARY 
TAD-44.1, ROOM 10402-B 
400 7TH STREET, S.W. 

WASHINGTON, D.C. 20590 
01CY ATTN R. LEWIS 
01CY ATTN R. DOHERTY 

INSTITUTE FOR TELECOM SCIENCES 
NATIONAL TELECOMMUNICATIONS & INFO AOMIN 
BOULDER, CO 80303 
01CY ATTN A. JEAN (UNCLASS ONLY) 

01CY ATTN W. UTLAUT 
01CY ATTN D. CROMBIE 
01CY ATTN L. BERRY 

NATIONAL OCEANIC & ATMOSPHERIC ADMIN 
ENVIRONMENTAL RESEARCH LABORATORIES 
DEPARTMENT OF COMMERCE 
BOULDER, CO 80302 
01CY ATTN R. GRUBB 
01CY ATTN AERONOMY LAB G. REID 
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DEPARTMENT OF DEFENSE CONTRACTORS 

AEROSPACE CORPORATION 
P.O. BOX 92957 
LOS ANGELES, CA 90009 
OICY ATTN I. GARFUNKEL 
01CY ATTN T. SALMI 
01CY ATTN V. JOSEPHSON 
OICY ATTN S. BOWER 
01CY ATTN N. STOCKWELL 
OICY ATTN D. OLSEN 
OICY ATTN J. CARTER 
OICY ATTN F. MORSE 
OICY ATTN SMFA FOR PWW 
OICY ATTN J. FENNEL 
OICY ATTN C. RICE 
OICY ATTN H. KOONS 

ANALYTICAL SYSTEMS ENGINEERING CORP 
5 OLD CONCORD ROAD 
BURLINGTON, MA 01803 
OICY ATTN RADIO SCIENCES 

BERKELEY RESEARCH ASSOCIATES, INC. 
P.O. BOX 983 
BERKELEY, CA 94701 
OICY ATTN J. WORKMAN 

BOEING COMPANY. THE 
P.O. BOX 3707 
SEATTLE, WA 98124 
OICY ATTN G. KEISTER 
OICY ATTN D. MURRAY 
OICY ATTN G. HALL 
OICY ATTN J. KENNEY 

CALIFORNIA AT SAN DIEGO. UNIV OF 
IPAPS, B-019 
LA JOLLA, CA 92093 
OICY ATTN HENRY G. BOOKER 
OICY ATTN E.C. WHIPPLE 

BROWN ENGINEERING COMPANY, INC. 
CUMMINGS RESEARCH PARK 
HUNTSVILLE, AL 35807 
OICY ATTN ROMEO A. OELIBERIS 

CHARLES STARK DRAPER LABORATORY, INC. 
555 TECHNOLOGY SQUARE 
CAMBRIDGE, MA 02139 
OICY ATTN D.B. COX 
OICY ATTN J.P. GILMORE 


COMPUTER SCIENCES CORPORATION 
6565 ARLINGTON BLVD 
FALLS CHURCH, VA 22046 
OICY ATTN H. BLANK 
OICY ATTN JOHN SPOOR 
OICY ATTN C. NAIL 

COMSAT LABORATORIES 
LINTHICUM ROAD 
CLARKSBURG, MD 20734 
OICY ATTN G. HYDE 

ELECTROSPACE SYSTEMS, INC. 

BOX 1359 

RICHARDSON, TX 75030 
OICY ATTN H. LOGSTON 
OICY ATTN SECURITY (PAUL PHILLIPS) 

ESL INCL. 

495 JAVA DRIVE 
SUNNYVALE, CA 94086 
OICY ATTN J. ROBERTS 
OICY ATTN JAMES MARSHALL 
OICY ATTN C.W. PRETTIE 


GENERAL ELECTRIC COMPANY 
SPACE DIVISION 
VALLEY FORGE SPACE CENTER 
GODDARD BLVD KING OF PRUSSIA 
P.O. BOX 8555 
PHILADELPHIA, PA 19101 
OICY ATTN M.H. BORTNER SPACE SCI LAB 

GENERAL ELECTRIC COMPANY 
P.O. BOX 1122 
SYRACUSE, NY 13201 
OICY ATTN F. REIBERT 

GENERAL ELECTRIC COMPANY 
TEMPO-CENTER FOR ADVANCED STUDIES 
816 STATE STREET (P.O. DRAWER QQ) 

SANTA BARBARA, CA 93102 
OICY ATTN DASIAC 
OICY ATTN DON CHANDLER 
OICY ATTN TOM BARRETT 
OICY ATTN TIM STEPHANS 
OICY ATTN WARREN S. KNAPP 
OICY ATTN WILLIAM MCNAMARA 
OICY ATTN B. GAMBILL 
OICY ATTN MACK STANTON 
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GENERAL ELECTRIC lECH SERVICES CO., INC. 

HMES 

COURT 

SYRACUSE, NY 13201 
01CY ATTN G. MILLMAN 

GENERAL RESEARCH CORPORATION 
SANTA BAR8ARA DIVISION 
P.O. BOX 6770 
SANTA 8ARBARA, CA 93111 
01CY ATTN JOE ISE JR 
01CY ATTN JOEL GARBARINO 

GEOPHYSICAL INSTITUTE 
UNIVERSITY OF ALASKA 
FAIRBANKS, AK 99701 
(ALL CLASS ATTN: SECURITY OFFICERS) 
01CY ATTN T.N. DAVIS (UNCL ONLY) 

OICY ATTN NEAL BROWN (UNCL ONLY) 

01CY ATTN TECHNICAL LIBRARY 
01CY ATTN T. HALLINAN 

ILLINOIS, UNIVERSITY OF 
ATTN: DAN McCLELLAN FOR K.C. YEH 
150 DAVENPORT HOUSE 
CHAMPAIGN, IL 61820 

INSTITUTE FOR DEFENSE ANALYSES 
400 ARMY-NAVY DRIVE 
ARLINGTON, VA 22202 
01CY ATTN J.M. AEIN 
01CY ATTN HANS WOLFHARO 
01CY ATTN JOEL BENGSTON 

HSS, INC. 

2 ALFRED CIRCLE 
BEOFORD, MA 01730 
01CY ATTN DONALD HANSEN 

INTL TEL & TELEGRAPH CORPORATION 
500 WASHINGTON AVENUE 
NUTLEY, NJ 07110 
01CY ATTN TECHNICAL LIBRARY 


JOHNS HOPKINS UNIVERSITY 
APPLIED PHYSICS LABORATORY 
JOHNS HOPKINS ROAD 
LAUREL, MD 20810 
01CY ATTN DOCUMENT LIBRARIAN 
01CY ATTN THOMAS POTEMRA 
OICY ATTN JOHN DASSOULAS 

LOCKHEED MISSILES 4 SPACE CO INC 
P.O. BOX 504 
SUNNYVALE, CA 94088 
OICY ATTN DEPT 60-12 
OICY ATTN D.R. CHURCHILL 

LOCKHEED MISSILES AND SPACE CO INC 
3251 HANOVER STREET 
PALO ALTO, CA 94304 
OICY ATTN MARTIN WALT DEPT 52-10 
OICY- ATTN RICHARO G. JOHNSON DEPT 52-12 
OICY ATTN W.L. IMHOF DEPT 52-12 
OICY ATTN 0. CAUFMAN 

KAMAN SCIENCES CORP 
P.O. BOX 7463 

COLORADO SPRINGS, CO 80933 
OICY ATTN T. MEAGHER 

LINKABIT CORP 
10453 ROSELLE 
SAN DIEGO, CA 92121 
OICY ATTN IRWIN JACOBS 
OICY ATTN I. ROTHMUELLER 

LOWELL RSCH FOUNDATION, UNIVERSITY OF 
450 AIKEN STREET 
LOWELL, MA 01854 
OICY ATTN K. BIBL 
OICY ATTN B. REINISCH 

M.I.T. LINCOLN LABORATORY 
P.O. BOX 73 
LEXINGTON, MA 02173 
OICY ATTN DAVID M. TOWLE 
OICY ATTN P. WALDRON 
OICY ATTN L. LOUGHLIN 
OICY ATTN D. CLARK 
OICY ATTN J. DAVIS 


JAYCOR 

1401 CAMINO DEL MAR 
OEL MAR, CA 92014 
OICY ATTN S.R. GOLDMAN 


MARTIN MARIETTA CORP 
ORLANDO DIVISION 
P.O. BOX 5837 
ORLANDO, FL 32805 
OICY ATTN R. HEFFNER 
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MCDONNELL DOUGLAS CORPORATION 
5301 BOLSA AVENUE 
HUNTINGTON BEACH, CA 92647 


OICY 

ATIN N. HARRIS 

OICY 

ATTN J. MOULE 

OICY 

ATTN GEORGE MROZ 

OICY 

ATTN W. OLSON 

OICY 

ATTN R.W. HALPRIN 

OICY 

ATTN TECHNICAL LIBRARY SERVICES 

MISSION 

RESEARCH CORPORATION 


735 STATE STREET 
SANTA BARBARA, CA 93101 
OICY ATTN P. FISCHER 
01CY ATTN W.F. CREVIER 
OICY ATTN STEVEN L. GUTSCHE 
01CY ATTN D. SAPPENFIELD 
01CY ATTN R. BOGUSCH 
01CY ATTN RALPH KILB 
01CY ATTN R. HENDRICK 
01CY ATTN DAVE SOWIE 
01CY ATTN F. FAJEN 
01CY ATTN M. SCHEIBE 
01CY ATTN CONRAD L. LONGMIRE 
01CY ATTN WARREN A. SCHLUETER 

MITRE CORPORATION, THE 
P.O. BOX 208 
BEDFORD, MA 01730 

01CY ATTN JOHN MORGANSTERN 
01CY ATTN G. HARDING 
01CY ATTN C.E. CALLAHAN 


PHOTOMETRICS, INC. 

44? MARRETT ROAD 
LEXINGTON, MA 02173 
01CY ATTN IRVING L. KOFSKY 

PHYSICAL DYNAMICS INC. 

P.O. BOX 3027 
BELLEVUE, WA 98009 
01CY ATTN E.J. FREMOUW 

PHYSICAL DYNAMICS INC. 

P.O. BOX 1069 
BERKELEY, CA 94701 
01CY ATTN A. THOMPSON 

R & D ASSOCIATES 
P.O. BOX 9695 
MARINA DEL REY, CA 90291 
01CY ATTN FORREST GILMORE 
01CY ATTN BRYAN GABBARD 
01CY ATTN WILLIAM B. WRIGHT JR. 
01CY ATTN WILLIAM J. KARZAS 
01CY ATTN ROBERT F. LELEVIER 
01CY ATTN H. ORY 
01CY ATTN C. MACDONALD 
01CY ATTN R. TURCO 

RAND CORPORATION, THE 
1700 MAIN STREET 
SANTA MONICA, CA 90406 
01CY ATTN CULLEN CRAIN 
01CY ATTN ED BEDROZIAN 


MITRE CORP 

WESTGATE RESEARCH PARK 
1820 DOLLY MADISON BLVD 
MCLEAN, VA 22101 
01CY ATTN W. HALL 
01CY ATTN W. FOSTER 

PACIFIC-SIERRA RESEARCH CORP 
1456 CLOVERFIELD BLVD. 

SANTA MONICA, CA 90404 
01CY ATTN E.C. FIELD JR 

PENNSYLVANIA STATE UNIVERSITY 
I0N1SPHERE RESEARCH LAB 
318 ELECTRICAL ENGINEERING EAST 
UNIVERSITY PARK, PA 16802 

(NO CLASSIFIED TO THIS ADDRESS) 

01CY ATTN IONOSPHERIC RESEARCH LAB 


RIVERSIDE RESEARCH INSTITUTE 
80 WEST END AVENUE 
NEW YORK, NY 10023 
OICY ATTN VINCE TRAPANI 

SCIENCE APPLICATION, INC. 

P.O. BOX 2351 
LAJOLLA, CA 92038 
OICY ATTN LEWIS M. LINSON 
OICY ATTN DANIEL A. HAMLIN 
OICY ATTN D. SACHS 
OICY ATTN E.A. STRAKER 
OICY ATTN CURTUS A. SMITH 
OICY ATTN JACK MCDOUGALL 

RAYTHEON CO. 

528 BOSTON POST ROAD 
SUDBURY, MA 01776 
OICY ATTN BARBARA ADAMS 
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SCIENCE APPLICATIONS, INC. 
HUNTSVILLE DIVISION 
2109 W. CLINTON AVENUE 
SUITE 700 

HUNTSVILLE, AL 35805 
01CY ATTN OALE H. OAVIS 

SCIENCE APPLICATIONS, INCORPORATED 
8400 WESTPARK DRIVE 
MCLEAN, VA 22101 
01CY ATTN J. COCKAYNE 

SCIENCE APPLICATIONS, INC. 

80 MISSION DRIVE 
PLEASANTON, CA 94566 
01CY ATTN SZ 

SRI INTERNATIONAL 
333 RAVENSWOOD AVENUE 
MENLO PARK, CA 94025 

01CY ATTN DONARD NEILSON 
01CY ATTN ALAN BURNS 
01CY ATTN G. SMITH 
01CY ATTN L.L. COBB 
01CY DAVID A. JOHNSON 
01CY ATTN WALTER G. CHESNUT 
01CY ATTN CHARLES L. RINO 
01CY ATTN WALTER JAYE 
01CY ATTN M. BARON 
01CY ATTN R. LIVINGSTON 
01CY ATTN RAY L. LEADABRAND 
01CY ATTN G. CARPENTER 
01CY ATTN G. PRICE 
01CY ATTN J. PETERSON 
01CY ATTN R. HAKE, JR. 

01CY ATTN V. GONZALES 
01CY ATTN 0. MCDANIEL 
01CY ATTN R. TSUNODA 

TECHNOLOGY INTERNATIONAL CORP 
75 WIGGINS AVENUE 
BEDFORD, MA 01730 
01CY ATTN W.P. BOQUIST 

UNIVERSITY OF TOKYO 
ISAS 

KOMABA, MEGURO-KU 
TOKYO, JAPAN 

01CY ATTN DR. K.I. OYAMA 


MAX-PLANCK-INSTITUT 
FUR PHYSIK UNO ASTROPHYSIK 
INSTITUT FUR EXTRATERRESTRICHE PHYSIK 
8046 GARCHING B. MUNCHEN, GERMANY 
01CY ATTN PROF. GERHARD HAERENDEL 


TRW DEFENSE 2 SPACE SYS GROUP 
ONE SPACE PARK 
REDONDO BEACH, CA 90278 
01CY ATTN R.K. PLE8UCH 
01CY ATTN S. ALTSCHULER 
01CY ATTN D. DEE 

VISIDYNE, INC. 

19 THIRD AVENUE 
NORTH WEST INDUSTRIAL PARK 
BURLINGTON, MA 01802 
01CY ATTN CHARLES HUMPHREY 
01CY ATTN J.W. CARPENTER 






IONOSPHERIC MODELING DISTRIBUTION i I ST 
UNCLASSIFIED ONLY 


PLEASE DISTRIBUTE ONE COPY (EXCEPT WHERE 
NOTED OTHERWISE) TO EACH OF THE FOLLOWING 
PEOPLE: 

ADVANCED RESEARCH PROJECTS AGENCY (ARPA) 
STRATEGIC TECHNOLOGY OFFICE 
ARLINGTON, VA 22217 

CAPT DONALD M. LEVINE 

NAVAL RESEARCH LABORATORY 
WASHINGTON, D.C. 20375 

DR. R. MEIER - CODE 4141 

DR. TIMOTHY COFFEY - CODE 4000 

DR. S. OSSAKOW - CODE 4780 

DR. J. GOODMAN - CODE 4180 
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TI1K STP/S3-4 SATELLITE EXPERIMENT: 
EQUATORIAL F-REGION IRREOULARI LIES 


INTRODUCTION 

Equatorial F-region irregularities have drawn much interest 
in recent years because of their undesirable effects on trans- 
ionospheric communications and their cause-effect relationship 
with fundamental plasma instability processes. In efforts to 
understand the causative mechanism(s), considerable advances have 
been made in areas of detailed ground-based radar observations 
and "in situ" measurements. The radar observations (e.g., 

Woodman and La Hoz, 1976; Tsunoda et al. 1979) have found that 
meter-size irregularities primarily populate the bottomside F- 
region in the early evening, at later times tend to rise up and 
break away from their lower altitude source regime and develop 
structures extending up to the 700-1000 km region. These structures 
have come to be called "plumes". 

On the macroscale (100's of meters to 10's of km) "in situ" 
nighttime equatorial measurements have revealed large biteouts in 
plasma density ranging up to three orders of magnitude (Hanson 
and Sanatani, 1973) and considered just as characteristic of 
spread-F as the much less intense (meter size) irregularities 
observed by radar. Later works on biteouts (Brinton et al. 1975; 
McClure et al. 1977) showed that ion composition inside and 
outside the holes can be vastly different. The molecular ions 
can be more abundant inside the holes than outside the holes, and 
the holes can vary from a few km to tens of km in the horizontal 
extent. 

In examining potential relationships between radar plumes 
and ionospheric depletions, Szuszczewicz (1973) suggested that 
equatorial holes and spread-F were the same phenomena with small 
scale irregularities imbedded within the large scale depletions. 

He argued that a chemical volume of ion density on the bottomside 
(containing the signature of bottomside species) could move 
upward through a stationary neutral atmosphere and appear as 
biteouts at higher altitudes with much smaller structures (down 
to the meter range) populating the density gradients which 
bounded the macroscale depletion. This model was in concert with 
the numerical results of Scannapieco and Ossakow (1976) and the 
drift measurements of McClure et al. (1977). 

To further study and definitively unfold the detailed 
relationships between large scale depletions, meter size 
irregularities and chemical transport processes, a coordinated 
investigation was conducted which involved simultaneous obser¬ 
vations by radar and "in situ" rocket-borne diagnostics 
(Szuszczewicz et al., 1980). The combined observations have 
shown that 

Manuscript submitted April 10, 1981. 
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(a) During conditions of well-developed equatorial spread-F 
the most intense "in situ" irregularities occurred on the bottom- 
side F-layer gradient. 

(b) Within a large scale topside F-layer depletion 
radar backscatter and "in situ" irregularity strengths 
maximized near the depletion's upper wall. 

(c) Ion composition within a topside depletion provided 
signatures of its bottomside source domain and estimates of 
average maximum vertical drift velocity. For long-lived 
depletions, it was found that molecular-ion signatures (N0 + 
and O 2 ) can be los£ while bottomside levels of N can be 
maintained when [0 ) ^ N g >> [NO ] + [0 2 ]; and finally, 

(d) Large scale fluctuations of 0 + accompanied by a 
near-constant level of NO and 0 2 on the bottomside F-layer 
gradient suggests that neutral atmospheric turbulence was 
not a major source for bottomside ionospheric plasma irregu¬ 
larities and the associated triggering of equatorial spread- 
F. 


To complement the vertical profile information provided 
by the rocket observations, we present and analyze a sample 
of "in situ" measurements conducted on the STP/S3-4 satellite 
carrying a pulsed plasma probe experiment. The probe experi¬ 
ment employed a self consistent test for measurement integrity, 
while determining electron density and temperature as well 
as density and mean ion mass fluctuations at 5-20 meter 
resolution. The S3-4 experiment has been discussed by 
Szuszczewicz et al. (1981). In this paper, we discuss some 
of the recent results with emphasis on general horizontal 
morphology, relationships to basic instability processes 
and associated scintillation effects. 

EXPERIMENTAL RESULTS 
Plasma Depletions 

The data reported here was made available by the NRL-747 
paired-pulsed-plasma-probes (P4) experiment (Szuszczewicz et 
al. 1981) on STP satellite S3-4 in a sun-synchronous orbit 
at lower F-region altitudes. The satellite crossed the 
nighttime equator at about 2230 LT when the the frequency of 
occurrence of spread-F was high. One probe was biased to 
respond to variations in plasma electron saturation currents 
(I « N ) and the other probe responded to the ion saturation 


2 











currents (I. “ N / \j M.) . Subject to the selection of a 
number of comman§able 1 modes of operation, either probe could 
be repetitively pulsed from its fixed-bias level using a 
special electronic procedure (Holmes and Szuszczewicz, 1975) 
to generate conventional Langmuir characteristics for full 
determinations of electron density N , temperature T and 
plasma potential V . The different modes of operation and 
experimental details are discussed by Szuszczewicz et al. 

(1981) . 

Figures 1 and 2 present samples of nighttime equatorial 
irregularity structures as measured by the currents collected 
by the electron (I ) and ion (I.) probes. The Figures show 
that the holes extlnd from a few kms to tens of kms, with 
depletion levels ranging from a factor of 3 (hole A) to a 
factor of 500 (hole F). A cursory analysis of Figures 1 and 
2 reveals three particularly interesting observations: 

(a) The density gradients on the opposing sides of each 
hole are different. Noting that the inclination of the 
satellite orbit is 96.4° (traveling in the east-to-west 
direction with time increasing left-to-right) it can be seen 
that density gradients in holes C-E (Fiq. 1) and F (Fiq. 2) 
are sharp on the eastern boundary and soft on each western 
counterpart. In the case of hole A, B and H (Fig. 1), the 
density gradients are not quite as sharp in comparison with 
holes C-F but the average density gradients are still softer 
on the western boundary. For the holes L (Fig. 1) and N' 
and L’ (Fig. 2), the difference in density gradients on the 
opposing sides does not constitute any specific behavior 
while in the hole C' (Fig. 2) the density gradient is sharp 

on the western boundary. These observations can be summarized 
by noting that the density gradients are different across 
the boundaries of each depletion, with a preference for the 
density gradients to be sharper on the eastern side. 

(b) Another interesting feature of the holes is their 
similarities in structural morphology. Holes A, B and H are 
similar in their horizontal profiles. A similarity feature 
can also be identified in the depletion L (Figure 1) and L' 
(Figure 2). Both of these depletions have a center point (M 
and M' respectively) around which the structures look similar. 
At the central point (M and M'), the plasma density approaches 
the background (undisturbed) level suggesting that one large 
hole is breaking up or alternatively two smaller ones are 
merging. 

(c) The top panel of Figures 1 and 2 displays relative 
density fluctuations, 6 I/I determined by variations about 
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PROBE CURRENT (AMP) 



LONG 34.2'E MLAT 3 6°S 

MLAT -10.4*S 

Figure 1 

Horizontal profile of ionospheric F-region plasma 
density as indicated by electron (I ) and ion (I.) 
saturation currents measured on rev e #2177. 6l/I 1 is 
the relative irregularity intensity as calculated 
by variations about linear detrends over sliding 
2.1 km intervals. 
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Figure 2 

Horizontal profile of ionospheric F-region plasma 
density as indicated by electron (I ) and ion (I.) 
saturation currents measured on rev e #2470. 61/I 1 

is the relative irregularity intensity as calculated 
by variations about linear detrends over sliding 
2.1 km intervals. 
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linear detrends over sliding 2.1 km intervals throughout the 
depletions. The fluctuations (6I/I) easily identities the 
degree of disturbance. The percent variation m the holes 
A, B and H (Figure 1) is less than 20% in most parts, while 
in other holes of Figures 1 and 2 the variation can be 
more than 50%. In addition, the variations can be different 
across the opposing boundaries of the holes. We discuss 
this more fully below. 


East-West Asymmetry 

Further illustrations of F-region irregularities appear 
in data revs 2122 and 2123 shown in Figure 3. Rev 2123 
shows four depletions (numbered 1 through 4) with the depletion 
level extending to more than two orders of magnitude. The 
rev 2122 shows three depletions (numbered 5 through 7) with 
depletions in density up to one order of magnitude. 

The discussion is facilitated by identifying certain 
features in rev 2123. First, there are clearly defined 
regions of undisturbed background ionosphere, marked alphabetically 
A through D; the smoothness of the relative density and the 
corresponding 0% fluctuations attest to their undisturbed 
nature. Focussing on depletions 3 and 4, we see that the 
irregularity intensities are 2 to 3 times larger on the 
western boundary than on eastern counterparts. This same 
relationship is true in depletions 1 and 2, but only after a 
qualification that suggests that 1 and 2 are halves of a 
larger depletion bounded by A and B. This is supported in 
part by the non-existence of a quiescient ionosphere between 
the two. When viewed from this perspective the western 
boundary is approximately twice as intense in irregularity 
intensity as the eastern boundary. In rev 2122, the depletions 
6 and 7 show the irregular intensity as more intense on the 
western boundary by a factor of 3-4. (On the basis of 
similar considerations applied to Figures 1 and 2, we note 
that the irregularity intensity is greater on the western 
boundaries of depletions D, E and H (Figure 1) and F (Figure 
2), while the eastern boundary in holes A-C (Figure 1) is 
more intense. Based on irregularity intensity (6I/I) 
observations, we conclude that there is a high probability 
of occurrence of more intense fluctuations across the western 
boundary of bottomside F-region depletions. 

The asymmetry in irregularity strengths and relationships 
to plasma instability mechanisms can be explored further 
through power spectral density analyses. We present in 
Figures 4 and 5 just such results for each of the boundaries 
in Figure 3 (IE and 1W refer respectively to the eastern and 
western boundary of depletion number 1). Though the experiment 
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provides the density fluctuation power spectra with a maximum 
Nyquist frequency of 400 Hz (19 m resolution along the orbit 
and 3-5 m resolution perpendicular to magnetic field) the 
data in this figure has been decimated by a factor of 3 in 
order to include greater spatial domains across the depletion 
walls. Power spectral analyses are presented across the 
boundaries of each of the seven depletions with spectral 
indices (n, in the equation P = P f n ) ranginq from 1.9 to 
2.5. More important however, is °the ratio of spectral strengths 
(entered as Pi hP q ( west)/P (east) in Figures 4 and 5) found to 
be 1.4 to 11.3 times more intense on the western boundaries. 

These intensity ratios extend down to a 15 meter wave¬ 
length perpendicular to the geomagnetic field. The spectral 
indices support the work of Keskinen et al. (1980) which 
predicts the same approximate range of values for horizontal 
irregularity structures perpendicular to B. The east- 
west asymmetry in the depletions is apparent from the 
linear detrend of data and the ratio of spectral strengths. 

We now show that the asymmetry bears on scintillation 
observations. To do this we note that scintillations depend 
upon AN(AN<*AI) rather than AN/N (as calculated by 

(Al) and the power spectral analysis P Q « -fA^) 2 ). For the 
pirpose of scintillation-effect calculations, we have determined 
|AI) r - m - s - the wavelength domain 80 m to 8 km by assuming 

that P r ' m ‘ s ' i n this domain is equal to (|Al| r ' m ‘ s '/I) 2 . 

Defining P 2 as 

p _ (AI) r m s (in the wavelength range 80rn-8km on the west wall) 

(AI) r m s tIie wa velength range 8 0m-8km on the east wall) 

we find the values of P 2 ranging from 1.1 to 7.7. These 
results suggest enhanced scintillation effects on the western 
wall of the depletions, in agreement with the observations 
of Livingston et al. (1980). Furthermore, the radar measure¬ 
ments (Tsunoda, 1979) also show that the bottomside back- 
scatter strength is often asymmetric in the east-west plane with 
stronger backscatter from the western wall of a plume. The 
combined observations support a model of E-W asymmetry which 
allows for a neutral-wind driven instability growth rate 
enhancement on the western side of a rising bottomside F- 
region depletion. The maximum growth rate occurs on the 
western wall of the rising depletion, where the electron 
density gradient is most closely aligned with the plasma 
drift velocity vector in the reference frame of the neutral 
wind (Tsunoda, 1979; Zalesak et al. 1980). 
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Figure U 


Power spectral analyses across the boundaries of the 
four depletions (1 through 4) of rev #2123 in Figure 3. 
The number in the upper right portion of each panel 
is the spectral index (n in the equation 

P * P Q k n ). The value P n in each right panel is the 
ratio of spectral strengths (P q ) of western to eastern 

boundary p = P (West)/P (East). The value P is the 
ratio of rTm.s. value of Spectral strength in the 
wavelength domain (80 m-8 Km) of western to eastern 

boundary (P 2 = p r . m . s . (West >/ p r . m . s .(East)). 
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Figure 5 

Power spectral analyses across the boundaries of the 
four depletions (5 through 7) of rev #2122 in Figure 3. 
The number in the upper right portion of each panel 
is the spectral index (n in the equation 

P = P k n )- The value P, in each right panel is the 
ratio 0 of spectral strengths (P q ) of western to eastern 

boundary P-^ = P Q (West)/P (East). The value P_ is the 
ratio of r.m.s. value or spectral strength in the 
wavelength domain (80 m-8 Km) of western to eastern 

boundary (P 0 = P (West)/P (East)). 

2 r.m.s. r.m.s. 








COMMENTS AND CONCLUSIONS 


The results of the high resolution SI-4 satellite 
experiment for equatorial F-region irregularities show large 
scale plasma depletions (1 km to 10's of km wide) with 
smaller scale irregularities superimposed (smallest detectable 
scale size ^ 10 meters perpendicular to B). The depletions 
(which may in fact represent the bottomside upwelling process 
that has come to be identified with the lower F-region 
manifestations of spread-F) show east-west asymmetry with 
the irregularity intensity and spectral strengths generally 
more intense on the western boundary. Associated calculations 
over density fluctuations in the range 80 meters to 8 km 
suggest that scintillation effects would be similarly more 
intense on the western boundary. We find these observations 
consistent with radar (Tsunoda, 1979) and scintillation 
measurements (Livingston et al. 1980) as well as the recent 
computational work of Zalesak et al. (1980). Furthermore, 
we find the horizontal power spectral indices to lie between 
1.9-2.5, supporting the role of Rayliegh-Taylor instability 
(Keskinen et al. 1980) and the recent theoretical and experi¬ 
mental comparison in the vertical plane (Keskinen et al. 

1981) . 
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nighttime ionosphere. Simultaneous measurements of 5M;/Mj -*• (Pjy| - A m f x m) and 6N e /N e -* (P N 
- AnfXn) have shown a general behavior tending to lower power (A m < A n ) and softer spectra (X^, 

< X„) in ion mass fluctuations when compared with fluctuations in total plasma density. Limited 
analyses of the two power spectral elements raises hopes for the differentiation between plasma 
mechanisms that can lead to similar indices in P^. 
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THE S3-4 IONOSPHERIC IRREGULARITIES SATELLITE EXPERIMENT: 

PROBE DETECTION OF MULTI-ION COMPONENT PLASMAS 
AND ASSOCIATED EFFECTS ON INSTABILITY PROCESSES 

INTRODUCTION 

In recent years interest in ionospheric research has 
turned from the definition of the zero-order ionosphere to 
an understanding of the more dynamic processes involving 
plasma instability mechanisms and their cause-effect relation¬ 
ships. While much progress has been made, more detailed 
experimental and theoretical work is needed to unfold the 
active first principles that govern the geoplasma environ¬ 
ment. For example, we have yet to define the electron 
density fluctuation power spectrum over the broad domain 
encompassed by equatorial spread—F (tens of kilometers to 

fractions of a meter). 

It is noted that density fluctuation power spectra 
represent an important approach to understanding the multi- 
stepped plasma processes in which large scale irregularities 

cascade to much smaller dimensions. In this regard the 
works of Dyson et al. (1974), along with the cataloging 

efforts of McClure and Hanson (1973) have Identified some 
characteristic spectral features down to 70 meters. From 
here the definition of spectral characteristics needs to be 
extended down to sizes less than 1 meter (e.g., Huba et al. 

1978 and associated references), and the same spectral 
characteristics must be studied within the context of positions 

relative to the F-layer peak (Ossakow et al, 1979) and any 
superimposed ionospheric depletions (McClure et al. 1977 and 
Szuszczewlcz, 1978). These measurements must be time correlated 
in order to study the development of multi-step instability 
processes. 
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Extended studies of fluctuation power spectra cannot 
stand alone in determining all the causal mechanisms. Existing 
equatorial data show that studies of F-region irregularities 
must determine the degree of balance between the chemistry 
and dynamics. To this end, the temperature and high resolution 
ion composition measurements must be made inside and outside 
of ionospheric holes, as well as across the sharp boundaries 
that are a characteristic feature of many of the depleted 
domains. The temperature measurement should also be done 
simultaneously with the determination of electron density 
power spectra in order to determine the role of electron 
energy in the naturally occurring instability processes 
(Ossakow, 1974) . 

In order to develop a more comprehensive experimental 
profile on ionospheric irregularity structures, the pulsed 
plasma probe technique (Szuszczewicz and Holmes, 1975) has 
been expanded to include measurements of electron density 
fluctuations and variations in mean ion mass. These measurements 
are conducted simultaneously with the standard determinations 
of absolute electron density (N^), temperature (T^) and 
plasma potential (V ). The capability for high spatial 
resolution measurements of mean ion mass opens up the possibility 
of exploring the coupling of electron density variations and 
ion composition irregularities by the comparison of electron 
density fluctuation power spectra P^(k) and a newly developed 
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parameter, the mean ion mass fluctuation power spectra, 

P^Ck). The new instrument (designated NRL-747) has been 
successfully flown on the lower F-region, polar orbiting, 
sun-synchronous (2230 LT equatorial crossing) STP/S3-4 
satellite launched in March 1978. The orbit allows a global 
study of N g , T^ and P^(k), an information set that will help 
catalogue similarities and differences between polar and 
equatorial irregularities and ultimately sort out first and 
second-order cause-effect relationships operating between 
plasma instabilities and ionospheric irregularities. The 
experiment provides the first effort to study the role of 
ion composition in the distribution of wave energy in the 
cascading process of large-to-small scale ionospheric plasma 
irregularities, a role considered important at the equator 
during the occurrence of spread-F conditions and potentially 
important at high latitudes under conditions of current- 
convective instability. 

Figure 1 schematically depicts the latitudinal differences 
in the irregular ionosphere at F-region altitudes, with one 
of the initial investigative targets in the S3-4 effort 
being the complete experimental definition of thermal electron 
distributions in the "biteouts" at nighttime equatorial 
latitudes under spread-F conditions. These "biteouts" are 
naturally occurring ionospheric holes which can be three 
decades deep in depletion and can have widths ranging from 
fractions of a kilometer to tens of kilometers (see e.g., 
McClure et al., 1973 and 1977; Szuszczewicz, 1978; and 








Szuszczewi.cz et al. , 1981). There can be major changes in 
ion chemistry which take place across the boundaries of the 
holes... changes which now appear to be trendable and strongly 
coupled to the mechanism(s) which generate the holes themselves 
(Szuszczewicz et al., 1980; Narclsi and Szuszczewicz, 1981), 
and cascade the irregularity distributions from 10's of 
kilometers to fractions of a meter (e.g., Keskinen et al.,1981; 
Kelly et al., 1981). 

Outside the latitudinal domain of equatorial spread-F, 
the regions of primary interest in the S3-4 study of F- 
region irregularities involve the main trough, auroral oval, 
and domains encompassing the ring current, polar wind and 
the cusp. The mid-latitude and dayside equatorial F-regions 
are generally very regular in structure and consequently of 
less interest. 

Within the framework of the relatively simple picture 
in Figure 1, the improved plasma diagnostic capabilities on 
the S3-4 satellite and associated analysis efforts will 
attempt to assemble a relatively comprehensive catalogue of 
polar and equatorial irregularities along with an improved 
understanding of their cause-effect relationships. Indeed, 
some initial efforts have already been published (Singh et 
al., 1981; Rodriguez et. al., 1981; and Szuszczewicz et al., 
1981). In the subsequent sections, the experimental technique 

will be described, associated data presented, and evidence 
will be provided which suggests that ion composition plays 

an important role in multi-stepped plasma instability 
processes active in the F-region ionosphere. 
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Fig. 1 — Schematic representation of F-region plasma density irregularities and associated geophysical domains 















THE EXPERIMENTAL TECHNIQUE 


The NRL-747 experiment employs a pair of pulsed plasma 
4 

probes (P is the designated acronym), each of which is 
capable of simultaneous measurement of electron density, 
temperature and density fluctuation power spectra regardless 
of the state of turbulence or the degree of irregularity in 
the ionospheric plasma medium. Together, the pair of probes 
also provide mean ion mass fluctuation measurements to a maximum 
Nyquist frequency of 200 Hz. The instrument is a Langmuir- 
type probe using a special electronic procedure for generating 
the current-voltage characteristic (Holmes and Szuszczewicz, 
1975; Szuszczewicz and Holmes, 1975). The pulse modulated 
approach reduces the distortion of the current-voltage 
characteristic that can result from surface contamination 
and allows millisecond tracking of density fluctuations that 
might occur during the time required to generate the current 
voltage characteristic. (The temporal resolution is limited 
only by telemetry.) 

Subject to the selection of 1 of 8 commandable modes of 
operation, each of the probes had applied to it some variation 
of the voltage function illustrated in Figure 2. The pulse 
modulated waveform, following the sawtooth envelope, provides 
the fundamental data set for a "conventional" Langmuir 
current-voltage characteristic and associated determinations 
of and T g (Chen, 1965) at a nominal 3 Hz rate. During 
the interpulse period, a fixed-voltage V_ is applied, with 
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associated current measurements providing a running measure 


of density fluctuations (assuming 5l_ « 5N ) and a time- 

B e 

dependent data sec for power spectral analysis with a Nyquist 

frequency of 400 Hz in the high data-rate mode. In addition 

to pulse-modulation, each of the probes could be operated in 

a simple fixed-bias mode and as a "conventional" continuous- 

sweep Langmuir probe, the latter having been introduced for 

systematic studies of hysteresis and contamination effects 

during satellite mission lifetimes (a minor mission objective). 

The probes were routinely operated with on one probe 

set for electron-saturation-current collection (defined as 

the E-probe with I„ = I (sat) = E) while the value of V„ on 
Be d 

the second probe was biased for ion saturation current 
collection (defined as the I probe with = I^sat) = I). 

The expressions for the currents collected by the two cylindrical 
probes take the forms 


, K.1 

E 3 I e (.«e) - N e V 


(Chen, 1965; for thick sheaths) and 



( 1 ) 



(Hoegy and Wharton, 1973; for probe axis perpendicular to 

the vehicle velocity vector in the ionospheric plasma rest 

frame). In the above equations A is the probe area, M ... 

P e(i) 

and N ... are the mass and density of the electron (ion) 
e ( i) 






population, T^^ is the associated temperature of an assumed 

Maxwellian distribution, e is the fundamental electron 

charge, k is Boltzmann's constant, w is the satellite velocity 
6 (i) 

and <f>p is the baseline voltage Vg applied to the E (I) 

probe and referenced to the plasma potential (<J> v / = V g - V # ). 
The square of the ratio I g (sat)/I.(sat) can be written 


, E _ l_ K I 1 + I e»p/« e I 

^7 IT j (M i w 2 /2.kT i + le^J/kTj 


with additional manipulation (assuming | e4>^ | >> kT g ) resulting in 


( I e (sat)> \ 

\ T^Tsat )- 1 


I (sat) 


$ 


2 i 

1 M . w < < e <b 
7 i P 


1 M. w2 >> ecj) 1 
7 1 P 


/ I e (sat) J 
y 1 1 (sat) / 


u 2 X -1 i M .» 2 % 


{ I 

2 °K / 


For laboratory and rocket-borne experiments equation (4a) 
would apply, whereas in the S3-4 satellite investigation 
equation (4c) applies. Equation (4c) is plotted in Figure 
3 for two sets of bias potentials, (|$p|, l^pl) “ (2V,1V) 
and (IV,2V). The results in Figure 3 show that with limitations 

q i 2 

placed on the ratio <t>^/<p^, variations in (I e /I^) can be 
taken as a direct measure of ion mass variations. 


9 







RATIO OF SATURATION CURRENTS (l e /l,) 



0 10 20 30 

ION MASS (AMU) 

Fig. 3 — Dependence of the saturation current ratio (I e /Ij) 2 on ion mass 









DATA ANALYSIS AND RESULTS 


Initial reduction routines for bulk processing and 
4 

plotting of P /S3-4 data begin with an orbit-by-orbit plot 
of relative electron density as measured by changes in ion- 
and electron-saturation currents. A representative sample 
of this data collected on orbit 2177 is shown in Figure 4 
where the abscissa coordinates are universal time (UT), 
altitude (ALT in km), latitude (LAT), longitude (LONG), 
magnetic latitude (MLAT), and L-shell value (L). The probes' 
magnetic aspect angle is also plotted in the figure. 

The left hand edge in Figure 4 corresponds to the 
satellite's ascending node (south-to-north) in the midnight 
hemisphere near the south magnetic pole (MLAT ■ -90°). 

With increasing time (UT) the satellite passes through the 
nighttime equator, the main trough, over the northern auroral 
oval and into the dayside ionosphere where vehicle solar cell 
voltage biases the entire vehicle such that both probes draw 
approximately equal ion saturation currents. 

The simultaneous measurements of electron- and ion- 
saturation currents, Ig(E) and I g (I) respectively, provide 
confidence that the observed irregularities involve plasma 
variations and not just secondary effects (e.g., aspect 
sensitivities and variations in spacecraft potential). 
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Figure 4, characteristic of data accumulated to date, 


shows that the nighttime F-region is orders of magnitude 
more irregular than its dayside counterpart. The irregularities 
in the nighttime southern high-latitude region (MLAT < - 40°) 
is considered characteristic of that domain, with the most 
intense structures generally showing electron density 
variations less than an order of magnitude. There are other 
characteristic features in the data within - 15° of the 
nighttime equator, where observations of large scale plasma 
depletions (1 to 3 orders of magnitude) occur with a frequency 
typical of equatorial spread-F. 

While data sets like chat shown in Figure 4 provide a 
global map of large scale ionospheric features , primary 
investigative objectives are directed at the relationships 
between the large scale features and much smaller scale 
irregularities (tens of meters and less) believed to result 

from multistepped plasma processes. To this end, the fundamental 

• 2 

data sets I (sat) and (I (sat)/I.(sat)) are Fast Fourier 
e el 

analyzed to determine their density fluctuation power 
spectra, P^(k) and P^(k) respectively, where 
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(5) 


and 




<VV 


6M, 


P H (k) 


( 6 ) 


The analytical relationship between 6N /N and 6M./M. can 

e e l l 

be simply established for a 2-component ion distribution of 

masses and densities (M , M„) and (N , N u ) respectively. Ue 

a p ap 

do this by use of the definitions 


M N + M Q N Q 
a a 6 6 


<V V 


(6a) 


N - N + N 
e a 6 


N a + N a 
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6N 


N 6 + N 8 
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+ Nr 
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(6b) 

(6c) 

(6d) 

(6e) 


and a straightforward manipulation to derive 


6M, 


dN 


f (a,6) 


(7a) 


where 


f(a,8) 



»>s 


B >8 


+ 1 


(7b) 
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It is appropriate to note that the experimental determination 

of mean-ion-mass fluctuations SM^(-*-P^), through variations 
2 

in [I (sat)/I. (sat)] , assumes the relative constancy of all 

e i 

potentials. (This includes the spacecraft potential as well 
as the potentials which each probe presents to the plasma.) 

The spacecraft potential can vary as a direct result of 
changes in local plasma density since the floating potential 
of a body is dependent upon the ratio of its radius to the 
Debye length. For large space vehicles however, floating 
potential variations caused by even substantial plasma 
density variations should be small (Szuszczewicz, 1972). 

Another possible source of potential variations Involve 
charging of contamination layers on the vehicle and/or on 

the probes (Szuszczewicz and Holmes, 1975). From the S3-4 

2 

data examined to date, variations in (I £ /I i ) associated 

with charging on contamination layers, appear to be a slowly 

varying function of time with no attendant effects on P w . 

Therefore, we have concluded that the spectral dependence of 

P M is Indeed representative of. variations in mean-ion-mass 
M 



To experimentally demonstrate P u (k) and the associated 

M 

relationship with P^(k) consider the high-resolution measure¬ 
ments (rev #2123) of the relative electron density across 

the nighttime equator (Figure 5). The peak electron density 

S—3 — 6 

is approximately 5 (10 ) cm at Ij ■ 3 (10 ) amp. The 

large scale depletions are seen to extend to 2 orders of 
magnitude with widths ranging from 50 to 170 km over a 600 km 
orbital segment. 
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UT 10:50:22 
ALT 241.1 km 
LAT-18.71 °S 
LONG 175.9°E 


TIME SECONDS 


10:58:42 
219 km 
11.9°N 
170.3°E 


Fig. & — An expanded view of relative electron density encountered during the nighttime equatorial 
crossing on S3-4 rev 2123. The relative electron density is presented by baseline electron saturation 
current. 


16 












In Figure 6 are presented the P„ and P„ results for a 

N M 

1 second interval located by point A in the density profile 
of Figure 5. Fitting the results to a power law behavior we 


f ind 


-2.9 

A f 
n 


-1.5 


(8a) 


and P ■ A f . (8b) 

M m 

By assuming that the time (frequency) domain spectral analysis 

in Figure 6 can be converted to wavelength through the 

* -2 9 

vehicle velocity (7.53 km/sec), the experiment shows f., 

N 

(<* k ^"^) from k 2rr/1 km 1 to k «* 2 tt/ 20 meters This 

is the first such satellite determination to wavelengths as 
short as 20 meters, with the earlier work of McClure and 
Hanson (1973) having defined some of the spectral features 
of equatorial spread-F down to 70 meters. (Conversion to 
the component of k perpendicular to the geomagnetic field 
extends the low wavelength end of Figure 6 down to k ■ 2ir/6 meter, 
the approximate value for the 0 + Larmor radius.) 

The spectral index for P N is approximately 15% steeper 
than previously reported values for conditions of bottomside 
spread-F (Keskinen et al., 1981) but well within the distribution 
of S3-4 spectral indices currently being accumulated and 
analyzed (Singh, currently unpublished) for conditions 
identified with the intermediate wavelength domain (k - 
2 tt/ 1 km" 1 to k - 2 tt/ 20 m" 1 ) . 


17 












METERS METERS 
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Fig. 6 Sample illustration of simultaneously determined density fluctuation power spectrum P^ and mean ion mass 

fluctuation power spectrum 











The P w « f observation is the first of its kind 

M 

and unique to the NRL-747/S3-4 experiment. Currently there 

are no computational guidelines on the expected behavior, 

but there is sufficient evidence in laboratory plasma studies 

to warrant such systematic considerations of ions and their 

role in the hierarchy of possible mechanisms covering the 

spectrum of observed ionospheric irregularities. The importance 

of ions is clear...even from the simple considerations of 

the Rayleigh-Taylor instability in which a difference in 

charged-particle drift velocities produces an electric field 

across a horizontal perturbation. These drift velocities 

are mass dependent (V^ « M^,(g x B)/B^) and vary directly as 

the mass of the i C ^ species. Similar mass discriminatory 

effects play an important role in ambipolar diffusion processes 

across gradients in plasma density. The process operates 

more rapidly on lighter ions and can result in "patches" of 

varying ion mass, with local variations in conductivity and 

electric fields, and ultimately an ion dependent Interaction 

in the process of energy dissipation in the large-to-small 

scale irregularity distribution. The P measurement has been 

M 

designed to test for just that type of interactive mode. 
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5Mi/Mi is a fairly complicated function of M^/Mg, 

N°/N° N^/ui and 6N /N itself (See Eq. 7b); and at this 
cl p ot p e e 

point we can only speculate on the many manifestations that 

P„ and P„ might take for the varied ionospheric conditions 
N M 

encountered in the S3-4 mission. For example, it has been 

suggested (J.R. Goldman, private communication) that differences 

in gradient scale lengths for N and li »’s would result in 

e a, p 

a more rapid fall off with increasing k for the quantity 

with an initially larger gradient scale length. This difference 

should be a direct observable through the P. T and P,. determination. 

Furthermore, there is the possibility that the simultaneous 

-2 

measurement of P (k) could help differentiate between a k 
H 

-2 

spectrum due to sharp edges and a k spectrum due to 
gradient-drift or say drift-dissipative waves. 

Currently, we have not analyzed sufficient data to 
exact any of the above analyses. We are at the beginning of 
a synoptic study designed to systematically unfold the 
interdependence of P^ and P^, and the contributing roles of 
ion mass distributions through M /M„ and N°/N°. However, 

up ctp 

the results in Figure 6 as well as in Figure 7 can be used 

to reflect some of our observations: 

(i) The absolute power in P N is generally much greater 

than that in P„. In Figure 6, for example, power law fits 
n 

yield A /A 86, reflecting a situation in which the 
7 n m 'v 

absolute density fluctuations 6N /N are two orders of 

e e 
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POWER SPECTRUM 


NRL 747 REV 2123 10:51:19 MODE B,E NRL 747 REV 2123 10:51:19 MODE B,(E/I) 



FREQUENCY (Hz) FREQUENCY (Hz) 


Fig. 7 — Four contiguous FFT’s of electron density fluctuations (left side A through D) and 
simultaneously determined fluctuations in mean ion mass (right side A through D) 











magnitude more intense than counterparts in 

(ii) More often than not the k-dependence of P„ is 

N 

maintained to shorter wavelength (higher frequencies) than P.. (k). 

M 

(iii) The relationship between.^the spectral indices of 

P N and P M is highly variable. Figure 7 is an illustration 

of a random selection of 1.2-second , interval spectra near 

the domain sampled in Figure 6. Panel A represents the power 

spectrum of the quiescent region before the depletion, while 

panels B through D represent contiguous power spectra within 

the depletion. In the quiescent region (Panel A), both 

density and mean ion mass fluctuations show noise spectra, 

while Inside the depletion well-defined k-dependencies in 

P^ and (less frequent) k-dependent spectra in P M are observed. 

While small in number, this sample of data is rather representative 

of results accumulated to date. More frequently than not P„ 

is equal to or softer than P N (e.g., panels B,C and D). In 

addition, there are observations of well-defined k-spectra 

in P N with none in P M (i.e., only noise, as in panel D). 

Panel D is interpreted as an unstable plasma environment 

with a single component ion population. This is in agreement 

with the predictions of Eq. (7b) with M q - Mg. 

Results like those outlined above are continuing with a 

view to a synoptic definition of the P , P interdependence 

n N 

and the associated roles of M^, gradient scale lengths and 

geoplasma conditions (N ,T and MLAT). 

e e 
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COMMENTS AND CONCLUSIONS 

In order to develop a more comprehensive experimental 
profile on ionospheric irregularities, the pulsed plasma 
probe technique has been expanded to include the measurement 
of mean ion mass and associated mass variations, 6M^/M^, on a 
spatial scale limited only by geometry and telemetry. Conducted 
simultaneously with the standard determinations of absolute 
density N e> temperature T , and density fluctuation power 
spectra dN^/N^ (k) , the expanded capability has opened 

up investigations which explore the coupling of electron and 
ion plasma wave energies which may influence the way in which 
large scale plasma perturbations cascade to much smaller 
dimensions. 

Although the relationship between 6M./M. and 6N /N 

xi e e 

is a complicated function of ratios in the ion populations, 
i.e., V M 6 , it may be systematically related 

to absolute densities, gradient scale lengths and time in 
cascading instability processes. 

The new technique has been successfully flown on an F- 
region, polar orbiting satellite with some of the early 
results including: 

(1) Extension of 5N measurement capability down to 
e 

dimensions of the order 5-20 meters, with associated power 
spectra of spectral index in the range 1.6 to 2.9. While 
short of the 1-3 meter limit desireable for comparison with 
Jicamarca or Altair radar coherent scatter results, the S3-4 


23 





experiment goes beyond the 70 meter limit of earlier satellite 
measurements and lends itself to additional arguments which 
support the association of the Rayleigh-Taylor instability 
occurrence of equatorial spread-F. 

(2) Measurements of 6M^ with resolution that allows 
detection of macroscale features (e.g., transitions from 0 + 
dominant to N0 + dominant environment) and microscale spectral 

analysis down to dimensions comparable to the 6N g capability. 

-X 

As expected, power law fits (P (k) ■ Af ) to P N and P^ have 
revealed considerable variations with X w generally (but not 
always) smaller than X^. Similarly, the absolute power A^ 
in ion mass fluctuations have been found to be considerably 
less (@ 2-100) than the power in P ^. A systematic study of 
the detailed relationships and dependence on zero-order 
plasma conditions is currently underway. 
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